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ABSTRACT 


The  oxidation  of  refractory  borides,  graphites  and  JT  composites, 
hypereutectic  carbide -graphite  composites,  refractory  metals,  coated 
rcfiaCtuxy  uiclala,  nxcLctl-uxide  composites,  and  iridium  coated  grapmtes 
in  air  has  been  studied  under  high  velocity  atmospheric  flight  conditions. 
Elucidation  of  the  relationship  between  hot  gas/cold  wall  (HG/CW)  and  cold 
gas/hot  wall  (CG/HW)  surface  effects  in  terms  of  heat  and  mass  transfer 
rates  at  high  temperatures  was  a  principal  goal. 

Published  arc  plasma  test  data  for  refractory  materials  taken  in  eight 
different  facilities  were  collected  and  examined  by  comparing  the  observed 
surface  temperature  with  calculated  radiation  equilibrium  values.  Wide 
variations  in  the  ratio  of  calculated  to  observed  temperature  were  encountered. 
Similar  calculations  performed  for  tests  conducted  in  the  present  program 
yielded  results  close  to  unity  especially  when  melting  is  encountered.  Larger 
ratios  (up  to  1.5)  were  noted  for  specific  materials  which  produce  silicon 
oxides,  implying  enhanced  resistance  to  energy  absorbtion.  Thus,  an  alter¬ 
native  method  of  presentation  which  compares  recession  rate  as  a  function  of 
heat  flux  and  enthalpy  for  the  candidate  materials  was  developed.  This  des- 
ci'iption  provides  a  means  for  comparing  performance  for  various  trajectories 
by  applying  a  flux/ enthalpy-altitude/velocity  translation  in  considering  candi¬ 
date  materials.  Comparison  of  the  trajectory  of  the  FDL-7MC  lifting  reentry 
vehicle  (Lift/Drag  ratio  between  2.  5  and  3.  0  and  a  3"  nose  radius)  eliminates 
all  of  the  candidate  materials  except  the  boride  composites.  These  composites 
have  survived  multicycle  exposures  totaling  20,000  seconds  under  conditions 
simulating  the  most  severe  portions  of  the  FDL-7MC  trajectory. 

Calculations  of  the  flux-enthalpy  boundaries  for  recession  rates  of  1 
mil/ sec  based  on  melting  of  the  solid  oxide  forming  materials  are  found  to 
compare  reasonably  with  observations.  The  model  employed  for  these  cal¬ 
culations  is  based  on  providing  the  latent  heat  required  for  melting  at  a  rate 
of  1  mil/ sec. 

Calculation  of  the  eurface  temperature  based  on  stream  and  material  proper¬ 
ties  is  presented  to  predict  internal  temperature  gradients  for  comparison  with 
the  "in-depth"  measurements.  Temperature  gradients  along  the  axis  of  a  right 
circular  cylinder  which  is  heated  from  one  end  in  an  arc  plasma  test  with  front 
face  and  side  radiation  losses  are  considered.  The  effects  of  radius,  length, 
thermophysical  properties  and  an  oxide  film  on  the  front  face  are  included. 
Measurement  of  temperature  gradients  through  oxide  films  formed  during  arc 
plasma  exposures  indicate  substantial  gradients  (1000°R  through  100  mils)  can 
exist.  Comparison  of  the  measurements  with  computed  results  yield  good 
agreement  in  view  of  the  simple  models  employed.  Explicit  models  are 
presented  for  computing  the  rate  of  graphite  recession  in  air  as  a  function  of 
density,  surface  temperature,  gas  velocity,  stagnation  pressure  and  sample 
radius.  The  results  are  compared  with  observations  covering  a  wide  range  of 
conditions . 


ThiB  abstract  is  subject  to  special  export  controls,  and  each  trans¬ 
mittal  to  foreign  governments  or  foreign  nationals  may  be  made  only 
with  prior  approval  of  the  Air  Force  Materials  Laboratory  (MAMC), 
Wright-Patterson  Air  Force  Base,  Ohio  45433. 


iii 


I 


9 


•7 


TABLE  OF  CONTENTS 


H  r\rt 


INTRODUCTION  AND  SUMMARY 


Introduction 


Summary 


THE  THERMAL  RESPONSE  OF  REFRACTORY  MATERIALS 
TO  HOT  GAS/COLD  WALL  EXPOSURES . 


Introduction 


Compilation  of  Experimental  Surface  Temperature 
and  Stream  Conditions  . . 


Correlation  of  Results 
Discussion  of  Results 


Results  Obtained  during  the  Present  Testing 
Program  . . . 


SIMPLIFIED  MODEL  FOR  CALCULATING  THE  OXIDATION 
BEHAVIOR  OF  GRAPHITE  UNDER  HIGH  VELOCITY  AIR 
FLOW  CONDITIONS  . 


Introduction  and  Summary 


B.  Derivation  of  the  Simplified  Model  for  Graphite 

Oxidation  16 

C.  Comparison  of  the  Scale- Gilbert  and  John-Schick 

Models  for  Graphite  Ablation .  20 

INFLUENCE  OF  FLUX- ENTHALPY  AND  ALTITUDE- 
VELOCITY  VARIABLES  ON  THE  RECESSION  OF  RE¬ 
FRACTORY  MATERIALS  22 

UTILIZATION  OF  THE  FLUX-ENTHALPY  VS.  ALTITUDE- 

VELOCITY  CORRELATION  TO  SCREEN  MATERIALS 

FOR  SPECIFIED  TRAJECTORIES .  25 

CALCULATION  OF  THE  FLUX- ENTHALPY  BOUNDARIES 
FOR  RECESSION  RATES  OF  1  MIL/SEC  VIA  MELTING  26 

CALCULATION  OF  TEMPERATURE  DISTRIBUTION  THROUGH 
A  COMPOSITE  CYLINDER  UNDER  STEADY  STATE 
CONDITIONS  ALLOWING  FOR  SIDE  LOSSES  VIA  RADIATION  28 


REFERENCES 


LIST  OF  TABLES 


Table  P-g 

1  Summary  of  Published  Flux -Temperature  Data  for 

Refractory  Materials  38 

2  Summary  of  Published  Flux- Temperature  Data  for 

Refractory  Materials  39 

3  Summary  of  Published  Flux- Temperature  Data  for 

Refractory  Materials  40 

4  Summary  of  Published  Flux- Temperature  Data  for 

Refractory  Materials  41 

5  Summary  of  Published  Flux-Temperature  Data  for 

Refractory  Materials  42 

6  Summary  of  Published  Flux- Temperature  Data  for 

Refractory  Materials  43 

7  Summary  of  Published  Flux-Temperature  Data  for 

Refractory  Materials  44 

8  Summary  of  Published  Flux- Temperature  Data  for 

Refractory  Materials  45 

9  Summary  of  Published  Flux-Temperature  Data  for 

Refractory  Materials  46 

10  Summary  of  Publii  *  Flux- Temperature  Data  for 

Refractory  Materials  47 

11  Summary  of  Published  Flux-Temperature  Data  for 

Refractory  Materials  48 

12  Summary  of  Published  Flux-Temperature  Data  for 

Refractory  Materials  49 

13  Summary  of  Published  Flux-Temperature  Data  for 

Refractory  Materials  50 

14  Summary  of  Published  Flux-Temperature  Data  for 

Refractory  Materials  51 

15  Summary  of  Published  Flux -Temperature  Data  for 

Refractory  Materials  52 

16  Summary  of  Published  Flux -Temperature  Data  for 

Refractory  Materials  53 


v 


Table 


Page 


i*7  T«»v!rv«ratiii*<»  anrf  Pressure  DeDendence  of  the  Enthalpy 
of  Air  (BTU/Lb) 

18  Averaged  Values  of  Total  Normal  Emittance  and  Ratio* 
of  Calculated  and  Observed  Temperatures  Derived  from 
Hot  Gas/Cold  Wall  Arc  Plasma  Tests 

19  Comparison  of  Calculated  and  Observed  Recession  Rates 
for  RVA(B5),  PT0178(B9),  POCO(B  10)  and  AT J  Graphites 
(Flat  Faced  Cylinders) 

20  Calculation  of  Flux-Enthalpy  Conditions  Corresponding 
to  the  Trajectory  for  the  FDL-7MC  Hypersonic  Reentry 
Vehicle  under  Maximum  Cross  Range  Conditions 


21  Summary  of  Data  used  in  Theoretical  Calculations  of 
Recession  Rates  as  a  Function  of  Flux-Enthalpy  Condi¬ 
tions 

22  Calculation  of  Temperature  Gradients  through  Arc 
Plasma  Test  Sample  ZrB2(A-3)-2MC 

23  Comparison  of  Observed  and  Calculated  Internal 
Temperature  Results 

24  Comparison  of  Observed  and  Calculated  Internal 
Temperature  Results 

25  Comparison  of  Observed  and  Calculated  Internal 
Temperature  Results 

26  n  Comparison  of  Observed  and  Calculated  Internal 

Temperature  Results 

27  Comparison  of  Observed  and  Calculated  Internal 
Temperature  Results 

28  Comparison  of  Observed  and  Calculated  Internal 
Temperature  Results 


61 

62 

63 

64 

65 

66 


67 


68 


LIST  OF  ILLUSTRATIONS 


Figure 

Page 

1 

Surface  Temperature-Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure.  Enthalpy  and  Surface  Emittance 

69 

2 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

70 

3 

Surface  Temperature-Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

71 

4 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

72 

5 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

73 

6 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

74 

7 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

75 

8 

Surface  Temperature-Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

76 

9 

Surface  Temperature-Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

77 

10 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

78 

11 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

79 

12 

Surface  Temperature-Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

80 

13 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

81 

14 

Surface  Temperature -He at  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

82 

15 

Surface  Temperature -Heat  Flux  Relations  as  a  Function 
of  Stagnation  Pressure,  Enthalpy  and  Surface  Emittance 

83 

16 

Ratio  of  Calculated  to  Observed  Surface  Temperatures 
Based  on  the  Cold  Wall  Heat  Transfer  Coefficient  and 
Radiation  Equilibrium  vs.  Cold  Wall  Heat  Flux 

84 

vii 


Figure 

Page 

17 

Ratio  of  Calculated  to  Observed  Surface  Temperatures 
Based  on  the  Cold  Wall  Heat  Transfer  Coefficient  and 
Radiation  Equilibrium  vs.  Stagnation  Enthalpy 

85 

in 

P.atio  of  Calculated  i,u  Observed  Surface  Temperatures 
Based  on  the  Cold  Wall  Heat  Transfer  Coefficient  and 
Radiation  Equilibrium  vs.  Stagnation  Pressure 

86 

19 

Ratio  of  Calculated  to  Observed  Surface  Temperatures 
Based  on  the  Cold  Wall  Heat  Transfer  Coefficient  and 
Radiation  Equilibrium  vs.  Mach  Number 

87 

20 

Ratio  of  Calculated  to  Observed  Surface  Temperatures 
Based  on  the  Cold  Wall  Heat  Transfer  Coefficient  and 
Radiation  Equilibrium  vs.  Observed  Temperature 

88 

21 

Ratio  of  Calculated  to  Observed  Surface  Temperatures 
Based  on  the  Fay-Riddell  Heat  Transfer  Coefficient 
and  Radiation  Equilibrium  vs.  Cold  Wall  Heat  Flux 

89 

22 

Comparison  of  Observed  and  Calculated  Recession 

Rates  for  RVA(B-5)  Graphite  Cylinders  at  3310°R 

90 

23 

Comparison  of  Observed  and  Calculated  Recession 

Rates  for  Graphite  Cones,  as  a  Function  of  Tempera* 
ture  and  Air  Velocity 

90 

24 

Calculated  Effects  of  Mach  Number  and  Body  Radius  on 
the  Oxidation  of  Graphite  at  Pe  a  1  atm 

91 

25 

Calculated  Effects  of  Mach  Number  and  Body  Radius  on 
the  Oxidation  of  Graphite  at  =0.1  at m 

92 

26 

Calculated  Effects  of  Mach  Number  and  Body  Radius  on 
the  Oxidation  of  Graphite  at  PQ  =  0.001  atm 

93 

27 

(Altitude/ Velocity)  vs.  Stagnation  Enthalpy /Flux)  Rela¬ 
tions  for  a  One  Inch  Sphere 

94 

28 

Recession  Rates  of  HfB 2  i(A-2),  ZrB2(A-3),  HfB»+20% 
SiC(A-4)(A-7)  and  H£B2  •f  35%SiC(A-9)  as  a  Function  of 
Heat  Flux  and  Stagnation  Enthalpy 

95 

29 

Recession  Rates  of  ZrB2+20%SiC(A-8)  and  ZrB2+l4% 
SiC+30%C(A-10)  as  a  Function  of  Heat  Flux  and  Stagna¬ 
tion  Enthalpy 

96 

30 

Recession  Rates  of  Boride  Z(A-5)  RVA(B-5),  PT0178 
(B-9)  and  POCO(B-IO)  as  a  Function  of  Heat  Flux  and 
Stagnation  Enthalpy 

97 

viii 


Figure 

Page 

31 

Recession  Rates  of  PG(B-6)  and  BPG(B-7)  as  a  Function 
of  Heat  Flux  and  Stagnation  Enthalpy 

98 

32 

Recession  Rates  of  HfC+ClC-1 11.  ZrC+C(C-12).  Si/RVC 
(B-8)  and  KT-SiC(E-14)  as  a  Function  of  Heat  Flux  and 
Stagnation  Enthalpy 

99 

33 

Recession  Rates  of  JTA(D-13),  JT0992(F-15)  and  JT0981 
(F-16)  as  a  Function  of  Heat  Flux  and  Stagnation  Enthalpy 

100 

34 

Recession  Rates  of  WSi2/W(G-18),  Sn-AJ/Ta-W(C-19), 
W+Zr+Cu(G-20)  and  W(Ag(G-21)  as  a  Function  of  Heat 

Flux  and  Stagnation  Enthalpy 

101 

35 

Recession  Rates  of  Si02+W(H-22)(H-23),  Hf-Ta-Mo(I-23) 
and  Ir/C(I-24)  as  a  Function  of  Heat  Flux  and  Stagnation 
Enthalpy 

102 

36 

Location  of  1  mil/ sec  Recession  Rate  Boundary  as  a 
Function  of  Heat  Flux  and  Enthalpy  for  Refractory 
Materials  at  Pe  =  1  atm. 

103 

37 

Recession  Rate  of  Hafnium  Diboride  as  a  Function  of 

Heat  Flux  and  Stagnation  Enthalpy 

104 

38 

Recession  Rate  of  Zirconium  Diboride  as  a  Function  of 
Heat  Flux  and  Stagnation  Enthalpy 

105 

39 

Recession  Rate  of  HfB2  +  SiG  Composites  as  a  Function 
of  Heat  Flux  and  Stagnation  Enthalpy 

106 

40 

Superposition  of  the  FDL-7MC  Trajectory  on  the  Flux- 
Enthalpy  Map  showing  Recession  Rates  for  Refractory 
Materials 

107 

41 

Computed  Boundaries  for  Melting  Recession  of  HfBjj  } 
(A-2) 

108 

42 

Computed  Boundaries  for  Melting  Recession  of  ZrBjj 
(A- 3) 

109 

43 

Computed  Boundaries  for  Melting  Recession  of  HfC+C 
(C-ll) 

110 

44 

Computed  Boundaries  for  Melting  Recession  of  ZrC+C 
(C-12) 

111 

45 

Variation  of  Tj  with  Tt  for  a  Number  of  Gs  Values 
According  to  Eq.  64  when  I/L  =  0  and  kp/kg  =0,10 

112 

ix 


Vvi  4  HHIu.VVVV.V'l' 


Figure 

Page 

46 

The  Effect  of  Interface  Temperature,  Tj,  on  the  Level 
of  the  Radiation  ContriK.^tier.,  Tp,  According  to  Eq.  66 

113 

47 

The  Dependence  of  the  Total  Radiation  Parameter 
(Fy+Fc)  on  ’f  j.  for  Various  V'aiues  of  Cg  when  I/L  = 

0. 10  and  kjr/kg  =  0. 10  Derived  from  Eqs.  66  and  69 

114 

48 

The  Effect  of  Cg  on  T*  and  the  Corresponding  Radiation 
Parameter  Fg  for  the  Case  where  I  =  0  (i.  e. ,  TpTf) 
Based  on  Eq,  74 

115 

X 


I, 


INTRODUCTION  AND  SUMMARY 


A.  Introduction 

The  response  of  refractory  materials  to  high  temperature 
oxidizing  conditions  imposed  by  furnace  heating  has  been  observed  to 
differ  markedly  from  the  behavior  in  arc  plasma  "reentry  simulators. " 

The  former  evaluations  are  normally  performed  for  lonjj;  times  at  fixed 
temperatures  and  slow  gas  flows  with  well -defined  solid/ gas -reactant/ 
product  chemistry.  The  latter  on  the  other  hand  are  usually  carried  out 
under  high  velocity  gas  flow  conditions  in  which  the  energy  flux  rather 
than  the  temperature  is  defined  and  significant  shear  forces  can  be  en¬ 
countered.  Consequently,  the  differences  in  philosophy,  observables, 
and  techniques  used  in  the  "material  centered"  regime  and  the  "environ¬ 
ment  centered,  reentry  simulation"  area  differ  so  significantly  as  to 
render  correlation  of  material  responses  at  high  and  low  speeds  difficult 
if  not  impossible  in  many  cases.  Under  these  circumstances,  expeditious 
utilization  of  the  vast  background  of  information  available  in  either  area 
for  optimum  matching  of  existing  material  systems  with  specific  missions 
or  prediction  and  synthesis  of  advanced  material  systems  to  meet  require¬ 
ments  of  projected  missions  is  sharply  curtailed. 

In  order  to  progress  toward  the  elimination  of  this  gap,  an 
integrated  study  of  the  response  of  refractory  materials  to  oxidation  in 
air  over  a  wide  range  of  time,  gas  velocity,  temperature  and  pressure 
has  been  designed  and  implemented.  This  interdisciplinary  study  spans 
the  heat  flux  and  boundary  layer -shear  spectrum  of  conditions  encountered 
during  high-velocity  atmospheric  flight  as  well  as  conditions  normally  em¬ 
ployed  in  conventional  materials  centered  investigations.  In  this  context, 
significant  efforts  have  been  directed  toward  elucidating  the  relationship 
between  hot  gas/cold  wall  HG/GW  and  cold  gas/hot  wall  CG/HW  surface 
effects  in  terms  of  heat  and  mass  transfer  rates  at  high  temperatures,  so 
that  full  utilization  of  both  types  of  experimental  data  can  be  made.  The 
elucidation  of  various  mass  transfer  reaction  regimes  have  been  studies  in 
gaseous  and  solid  oxide  formation. 

The  principal  goal  of  this  study  is  the  coupling  of  the  material- 
centered  and  environment- centered  philosophies  in  order  to  gain  abetter  in¬ 
sight  into  systems  behavior  under  high-speed  atmospheric  flight  conditions. 
This  coupling  function  has  been  provided  by  an  interdisciplinary  panel  com¬ 
posed  of  scientists  representing  the  component  philosophies.  The  coupling 
framework  consists  of  an  intimate  mixture  of  theoretical  and  experimental 
studies  specifically  designed  to  overlap  temperature/energy  and  pressure/ 
velocity  conditions.  This  overlap  has  provided  a  means  for  the  evaluation 
of  test  techniques  and  the  performance  of  specific  materials  systems  under 
a  wide  range  of  flight  conditions.  In  addition,  it  provides  abase  for  develop¬ 
ing  an  integrated  theory  or  modus  operandi  capable  of  translating  reentry 
systems  requirements  such  as  velocity,  altitude,  configuration,  and  life 
time  into  requisite  materials  properties  as  vaporization  rates,  oxidation 
kinetics,  density,  etc.,  over  a  wide  range  of  conditions. 
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The  correlation  of  heat  flux,  stagnation  enthalpy,  Mach  No.  , 
stagnation  pressure,  and  specimen  geometry  with  surface  temperature 
through  the  utilization  of  thermodynamic,  thermal  and  radiational  proper¬ 
ties  of  the  material  and  environmental  systems  used  in  this  study  was  of 
prime  importance  in  defining  the  conditions  for  overlap  between  materials- 
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Significant  practical  as  well  as  fundamental  progress  along  the 
above  uiemiuued  lines  necessitated  evaluation  of  refractory  material  systems 
which  exhibit  varying  gradations  of  stability  above  2700°F.  Emphasis  was 
placed  on  candidates  for  3400  F  to  6000°F  exploitation.  Thus,  borides, 
carbides,  boride- graphite  composites  (JTA),  JT  composites,  carbide- 
graphite  composites,  pyrolytic  and  bulk  graphite,  PT  graphite,  coated 
refractory  metals/alloys,  oxide-metal  composites ,  oxidation- resistant 
refractory  metal  alloys,  and  coated  graphites  were  considered.  Similarly, 
a  range  of  test  facilities  and  techniques  including  oxygen  pickup  measure¬ 
ments,  cold  sample  hot  gas,  and  hot  sample  cold  gas  devices  at  low  velocities 
as  well  as  different  arc  plasma  facilities  capable  of  covering  the  50-2500 
BTU/ft^sec  flux  range  under  conditions  equivalent  to  speeds  up  to  Mach  12 
at  altitudes  up  to  200,000  ft  were  employed.  Stagnation  pressures  covered 
the  range  between  0.001  and  10  atmospheres.  .  Splash  and  pipe  tests  were 
performed  in  order  to  evaluate  the  effects  of  aerodynamic  shear.  Based  on 
the  present  results,  this  range  of  heat  flux  and  stagnation  enthalpy  produced 
surface  temperatures  between  2000°F  and  6500°F. 

B.  Summary 

The  present  report,  which  is  the  seventh  in  a  series  (_l-6)*,  deals 
with  theoretical  methods  for  correlating  the  performance  of  the  candidate 
refractory  materials  with  stream  conditions.  In  many  respects,  this  correla¬ 
tion  constitutes  the  essence  of  the  entire  study.  Thus,  given  stream  character¬ 
istics  such  as  stagnation  enthalpy,  stagnation  pressure  and  cold  wall  heat  flux, 
is  it  possible  to  provide  a  means  for  predicting  the  response  of  a  candidate 
material  upon  insertion  into  this  Stream?  The  development  of  such  methods 
would  be  of  twofold  value  to  the  current  program.  Operationally  it  would 
provide  a  means  for  checking  the  internal  consistency  of  arc  plasma  tests. 

Thus,  measurements  of  arc  plasma  stream  conditions  could  be  employed  to 
compute  the  surface  temperature  achieved  by  candidate  materials  during 
exposure  for  comparison  with  observed  surface  temperature.  In  addition,  since 
the  behavior  of  these  materials  is  strongly  dependent  on  the  surface 
temperature  (or  temperature  regime)  the  foregoing  correlation  could  provide 
a  means  of  describing  the  materials  performance  in  terms  of  flight  character¬ 
istics.  This  translation  could  be  affected  by  employing  the  relations  between 
altitude,  velocity  and  body  radius  on  the  one  hand  and  enthalpy,  pressure  and 
heat  flux  on  the  other.  In  this  way,  a  logical  method  for  comparing  the  require¬ 
ments  of  specific  flight  trajectories  with  the  capabilities  of  refractory  materials 
could  be  developed. 

Accordingly,  activities  aimed  at  generating  such  methods  were 
carried  out  during  the  course  of  the  program.  At  the  outset,  a  literature 
survey  of  published  arc  plasma  test  data  for  refractory  materials  was  per¬ 
formed.  Data  taken  in  eight  different  facilities  were  collected  and  examined 
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by  comparing  the  observed  surface  temperature  with  values  calculated  on 
the  basis  of  radiation  equilibrium.  Rather  wide  variations  in  the  ratio  of 
calculated  to  observed  temperature,  T(GALC)/T(OBS)  were  encountered. 

In  most  'nstances  this  ratio  was  greater  than  unity  and  in  some  instances 
(exposures  performed  in  a  given  facility)  the  ratio  T(GALC)/T(OBS)  re¬ 
sulted  in  values  near  2,0.  Similar  icjlaticns  performed  tor  all  of 

ihe  arc  plasma  tests  conducted  in  the  present  program  (6).  In  these  cases 
(nearly  800  in  total)  the  ratio  T(CALC)/T(OBS)  was  reasonably  close  to 
unity.  In  particular  the  present  results  indicate  that  values  near  unity 
are  observed  when  melting  is  encountered.  Larger  ratios  (up  to  1.5)  were 
noted  for  specific  materials  which  produce  silicon  oxides  such  as  H£B£+ 
SiC(A-4),  KT-SiC(E-14),  and  WSiz/W  (G-18)  and  for  Sn-Al/Ta-10W(G-19).  At 
low  temperatures  (i.  e. ,  3000^-3500^)  these  materials  exhibit  T(CALC)/ 
T(OBS)  ratios  near  1.5  providing  that  melting  does  not  occur.  The  occur¬ 
rence  of  ratios  which  are  larger  than  unity  implies  enhanced  temperature 
capability  due  to  resistance  to  energy  ahaorbtion  by  the  material.  Although 
the  origin  of  this  resistance  is  not  clear  at  present,  it  is  probably  due  to 
blocking  effects  caused  by  evolution  of  gaseous  oxides.  These  observations 
suggest  a  method  of  ranking  the  behavior  of  the  refractory  materials  which 
differs  from  the  customary  recession  vs.  temperature  curves.  Thus,  an 
alternative  method  of  presentation  which  compares  recession  rate  as  a 
function  of  heat  flux  and  enthalpy  for  the  candidate  materials  was  developed. 
This  method  does  not  require  a  knowledge  of  the  spectral  or  the  normal  emit- 
taace  and  integrates  the  blocking  effects  characteristic  of  each  material. 

In  the  course  of  the  present  study,  the  oxidation  of  graphites  in 
air  has  been  investigated  experimentally  over  a  range  of  conditions  (4-6) 
between  2500°R  and  6500°R,  at  velocities  between  1  ft/ sec  and  Mach  3.2. 

The  succeeding  volume  of  this  series  presents  a  complete  discussion  of  the 
surface  reaction  problem  encountered  in  the  oxidation  of  graphite.  This 
discussion  considers  the  coupling  of  mass  transport  through  the  boundary 
layer  with  reactions  at  the  surface  in  detail.  By  contrast,  the  discussion 
presented  here  employs  simplified  models  which  provide  an  explicit  means 
for  computing  the  rate  of  graphite  recession  in  air  as  a  function  of  density, 
surface  temperature,  gas  velocity,  stagnation  pressure  and  sample  radius. 
The  results  are  compared  with  observations  covering  a  range  of  (tensity 
between  80  and  115  lbs/ft^,  temperatures  between  2500°  and  6500  R, 
velocities  between  1  ft/ sec  and  Mach  8.  0,  stagnation  pressures  between  0.  007 
and  1.0  atm  and  nose  radii  between  0. 005  and  0.  07  ft.  The  current  descrip¬ 
tion  is  based  on  the  product  of  Arrhenius'  term  and  an  oxygen  partial  pres¬ 
sure  term.  The  former  consists  of  a  pre-exponential  of  0.74  lbs /ft^ sec  and 
an  activation  energy  of  10,730  cal/mole.  The  oxygen  partial  pressure  term 
has  an  exponent  of  0. 333  and  is  modified  by  an  explicit  correction  factor 
which  relates  the  oxygen  concentration  at  the  reacting  surface  to  the  oxygen 
concentration  at  the  edge  of  the  boundary  layer.  This  correction  factor  is 
specified  in  terms  of  Mach  No. ,  body  radius  and  pressure. 

A  method  for  describing  the  response  of  refractory  materials  to 
the  enthalpy  and  heat  flux  characteristics  of  the  stream  has  been  developed. 
This  description  provides  ameaas  for  comparing  material  performance  at 
the  stagnation  point  for  various  trajectories.  The  comparison  can  be  made 
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by  translating  the  flux- enthalpy  description  into  altitude- velocity  character¬ 
istics  based  on  established  relations  between  stagnation  pressure,  altitude 
and  velocity.  The  Fay-Riddell  relation  is  employed  to  specify  heat  flux  in 
terms  of  altitude,  velocity  and  body  radius.  The  material  ranking  afforded 
by  this  description  shows  that  HfB2+SiC{A-4)  possesses  the  widest  range  of 
applicability  of  all  the  candidate  materials  investigated  in  the  present  study. 

In  addition,  the  different  mcdcc  of  behavior.  cjJjJ.uii.eci  by  ablators  (such  as 
graphite  and  tungsten)  and  solid  oxide  formers  are  clearly  displayed. 

The  applicability  of  the  fiux/enthalpy-altitude/velocity  descrip¬ 
tion  in  considering  candidate  materials  has  been  illustrated  by  comparing 
the  trajectory  of  the  FDL-7MC  lifting  reentry  vehicle  with  the  behavior  of 
candidate  refractory  materials.  This  vehicle  is  designed  for  a  Lift/jJrag 
ratio  between  2.5  and  3.0.  The  conditions  imposed  by  this  trajectory  for 
the  case  of  a  3"  nose  radius  eliminate  all  of  the  candidate  materials  except 
the  boride  composites.  Theue  composites  have  survived  multicycle  expo¬ 
sures  totaling  20, 00 0  seconds  under  conditions  simulating  the  most  severe 
portions  of  the  FDL-7MC  trajectory  (6). 

Calculations  of  the  flux- enthalpy  boundaries  for  recession  rates 
of  1  mil/sec  based  on  melting  of  the  solid  oxide  forming  materials  are  found 
to  compare  reasonably  with  observations.  The  model  employed  for  these 
calculations  is  based  on  providing  the  latent  heat  required  for  melting  at  a 
rate  of  1  mil/sec. 

Measurement  of  temperature  gradients  which  exist  through 
oxide  films  formed  during  arc  plasma  exposures  indicate  substantial 
gradients  (1000°R  through  100  mils)  can  exist  (6). 

A  first  order  calculation  of  the  surface  temperature  as  a  function 
of  stream  conditions  and  material  properties  is  presented  in  order  to  pro¬ 
vide  a  means  for  predicting  internal  temperature  gradients  for  comparison 
with  the  "in  depth"  measurements.  This  calculation  considers  temperature 
gradients  along  the  axis  of  a  right  circular  cylinder  which  is  heated  from  one 
end  in  an  arc  plasma  test.  Front  face  and  side  radiation  losses  are  con¬ 
sidered  in  describing  the  effects  of  radius,  length  and  thermophysical  pro¬ 
perties  on  the  surface  temperature  and  internal  gradients.  In  addition,  the 
effect  of  an  oxide  film  on  the  front  face  was  included.  The  calculations  in¬ 
dicate  that  small  gradients  occur  when  the  surface  temperature  is  low,  or 
when  the  cylinder  length  and/or  the  length/radius  ratio  is  small.  Large  values 
of  the  thermal  conductivity  of  the  cylinder  material  also  leads  to  small  gra¬ 
dients.  Large  values  of  surface  temperature,  length/radius  ratio  and  small 
values  of  the  thermal  conductivity  of  the  oxide  and  base  material  re  suit  in  large 
gradients. 


The  model  has  been  applied  to  calculation  of  temperature  gra¬ 
dients  for  comparison  with  the  experimental  results  obtained  in  sixty- five 
arc  plasma  tests  on  a  variety  of  refractory  materials.  Ablators  and  oxide 
forming  materials  covering  a  wide  range  of  thermopbysical  and  oxidation 
characteristics,  Buch  as  ZrB2+SiC(A-8),  ZrB2(A-3),  HfB2+SiC(A-7),  RVA 
(B-5),  ZrB2+SiC+C(A-10),  WSi2/W(G-18)  and  Hf-Ta-Mo(I-23)  were  included. 
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Observables  consisted  of  the  measured  front  face  temperature  Tr,  the  ob¬ 
served  temperature  at  a  distance,  d  mils  from  the  front  face,  T(d),  the 
cold  wall  heat  flux,  q,  the  stagnation  enthalpy,  ie,  and  the  stagnation  pres¬ 
sure  Pe.  Additional  input  consisted  of  sample  radius,  R,  length,  L,  and 
oxide  coating  thickness,  1.  The  latter  was  equated  to  the  conversion  depth 
for  the  oxide  formers  (6).  For  WS^/W,  I  was  equated  to  the  WSi2  coat¬ 
ing  thickness  with  1=0  for  RVA(B-5)  graphite  which  ablates  without  coating 
formation.  Suitable  values  of  the  emittance,  e  g,  and  the  thermal  conductiv¬ 
ities  of  the  coating,  kp,  and  the  substrate,  kg,  Vre re  also  employed. 

The  computed  results  are  displayed  in  terms  of  the  ratio  of 
calculated  front  face  temperature  to  observed  front  face  temperature 
Tf<CALC)/Tf(OBS)  and  the  ratio  of  computed  in-depth  temperature  T^{CALC) 
to  computed  front  face  temperature  Tf(CAIiC).  The  latter  is  compared  with 
the  ratio  of  observed  in-depth  temperature  T^(OBS)  to  observed  front  face 
temperature  Tf(OBS).  If  agreement  between  calculated  and  observed 
temperatures  results,  TdCALC )/ lAOBS)  would  equal  unity  and  the  ratios 
Td(CAljC)/Tf(CALC)  and  Td(OBS)/Tf(OBS)  would  coincide.  Relatively  good 
agreement  was  encountered  as  regards  the  latter  comparison  in  view  of  the 
simple  model  employed  to  describe  the  complex  tests.  In  addition, 
Tf(CALC)/Tf(OBS)  ratios  were  computed  near  unity  for  many  of  these  tests. 
However,  in  line  with  the  behavior  noted  above  some  systematic  deviations 
were  observed. 

The  largest  of  these  occurred  at  low  surface  temperatures 
(i.e. ,  Tf  <  3300°R)  for  the  materials  which  form  SiC>2  as  an  oxidation 
product.  Thus,  in  cases  where  samples  of  HfB2+SiC(A-7),  ZrB2+SiC(A-8), 
ZrB2+SiC+C(A- 10)  or  WSi2/w(G-18)  were  exposed  with  shrouds  or  as  large 
diameter  hemispheres,  Tf(CALC)  is  considerably  larger  than  Tf(GBS).  How¬ 
ever,  this  difference  is  smaller  than  obtained  when  T|  is  computed  on  the  basis 
of  front  face  radiation  equilibrium  alone.  The  cause  of  this  behavior  is 
presently  unknown. 
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THE  THERMAL,  RESPONSE  OF  REFRACTORY  MATERIALS  TO 
fft5T"GA57 COLD  WALL  EXPOSURES 


A.  Introduction 


Prior  to  launching  the  extensive  HG/ CW  testing  program  con¬ 
ducted  under  the  presen:  investigation,  published  data  on  the  surface 
temperature  of  borides,  graphites,  graphite  composites,  silicon  carbide, 
boron  nitride,  tungsten  alloys  and  composites,  refractory  metal-oxide 
composites,coated  refractory  metals  and  iridium  coated  graphite  ex¬ 
posed  in  arc  plasma  and  Wave  Superheater  tests  were  collected.  These 
data  covered  exposure  conditions  over  a  range  of  stagnation  pressures 
betweenO.  002  and  70  atmospheres  while  Mach  Numbers,  stagnation 
enthalpy  and  heat  flux  levels  ran  from  0. 2  to  8. 5,  1400  to  18, 000  BTU/lb 
and  20  to  4080  BTU/ft^sec,  respectively.  The  specimen  configurations 
tested  included  flat  face  cylinders  and  hemispherical  caps  with  diameters 
between  0.25  and  3.00  inches.  The  facilities  at  which  these  exposures 
were  performed  included  Avco/SSD,  Cornell  Aeronautical  Laboratory, 
General  Electric  Space  Science  Center,  Plasmadyne  Corporation, 
Cinncinati  Testing  Laboratory,  North  American  Aviation  Center, 
Grumman  Aircraft  Engineering  Corporation  and  Aerospace  Corporation. 
The  results  were  compared  with  computed  surface  temperatures  based 
on  radiation  equilibrium  in  order  to  estimate  the  current  level  of  reli¬ 
ability  of  surface  temperature  predictions  from  stream  conditions.  A 
similar  comparison  of  the  re  suit  8  obtained  in  the  present  program  (6) 
has  been  provided  and  is  compared  with  the  earlier  results. 

B.  Compilation  of  Experimental  Surface  Temperatures  and 

Stream  Conditions  “ 

Since  the  present  study  is  concerned  with  a  comparison  of 
CC/HW  and  HG/CW  oxidation,  correlation  of  the  surface  temperature 
of  models  (specimens)  exposed  under  HG/CW  conditions  with  specific 
stream  characteristics  is  of  paramount  importance.  In  order  to  gain 
some  insight  into  the  relationship  between  the  surface  temperature  and 
stagnation  enthalpy,  heat  flux  and  pressure  under  subsonic  and  super¬ 
sonic  flow  conditions,  a  review  of  available  literature  has  been  made. 

The  results  are  contained  in  Tables  1-16  which  identifies  the  material, 
Mach  Number,  stagnation  pressure  (Pe),  enthalpy  (L),  cold  wall  heat 
flux  (q^)  and  model  configuration  (D) .  As  indicated^  the  diameter  of 
cylindrical  samples  are  designated  by  asterisks.  In  addition,  the  ob¬ 
served  surface  temperature,  measured  optically,  is  noted.  Since  the 
assignment  of  ‘'observed  temperature"  is  not  performed  uniformly  in 
Tables  1-16,  it  is  worthwhile  to  report  the  methods  employed  for  each 
set  of  measurements.  No  attempt  will  be  made  at  present  to  correct 
these  "observed  temperatures". 

The  ZrB2>  HfB£  and  H£B2~SiC  exposures  reported  under 
Reference  (7)  in  Tables  1  and  2  were  performed  at  Avco/SSD.  Bright¬ 
ness  temperatures  were  measured  at  X  =0.  65  microns  and  converted 
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to  "observed  temperature"  by  using  an  emittance  value  of  0.60.  Exposure 
times  of  five  to  thirty  minutes  were  employed.  The  high  pressure  ex¬ 
posures  of  ZrB 2,  Boride  Z,  ZrC^,  etc. ,  shown  under  Reference  (8) 
in  Tables  2-8  were  performed  in  the  Cornell  Wave  Superheater.  The 
former  runs,  as  well  as  those  in  Tables  4-6  and  8  designated  by  Refer¬ 
ence  (11)  (which  =* I?"  carried  out  iu  the  Wave  Superbeater),  were 

given!?  second  exposure  times.  In  these  cases,  the  brightness  tempera¬ 
ture  at  V  =  2.1  +  0.5  microns  is  reported.  These  temperatures  represent 
mean  values  observed  during  exposure. 

The  graphite  an.?  cugsten  exposures  in  Tables  7,  9,  12  and  13 
denoted  by  Reference  (10)  were  performed  at  Avco/feSD.  These  tests  were 
of  60  to  120  seconds  duration.  Total  emittance  values  were  measured  with 
an  Eppley  thermopile.  In  addition,  brightness  temperatures  were  measured 
at  h  a  0.65  microns  and  converted  to  surface  temperature  by  assuming 
that  the  emittance  at  X  =  0.65  microns  is  equal  to  the  total  emittance. 

The  results  reported  in  Table  2  for  tantalum  (9)  were  obtained 
during  50-150  second  exposures  using  an  optical  pyrometer  at  X.  =2.1  + 

0.4  microns.  Brightness  temperatures  weire  converted  to  observed  tem¬ 
peratures  by  using  an  emittance  of  0.45.  The  latter  value  was  obtained 
from  the  observation  of  the  melting  point  of  TajOg  at  3730°R.  The  HfC- 
C  (20)  exposure  time  was  6  seconds.  The  temperature  level  of  this  CAL. 
Wave  Superheater  test  was  established  by  observation  of  incipient  melting 
of  HfOz  (Table  4). 

The  ATJ  Graphite  exposures  denoted  by  Reference  (12)  were 
performed  in  the  General  Electric  tandem  Gerdien  and  free  jet  facilities 
(Table  10).  The  RVA  graphite  and  graphite  composites  designated  by 
Reference  (13),  Tables  11  and  12  were  similarly  tested.  Exposure  time 
was  60  to  1000  seconds  and  surface  temperatures  were  measured  with  a 
two  color  pyrometer. 

Exposure  times  of  100-1200  seconds  were  employed  testing 
the  JTA  and  Zr02  materials  shown  in  Tables  13  and  14  designated  by 
Reference  (16).  The  reported  surface  temperatures  are  optical  brightness 
temperature  measured  at  X  =  0.65p.  The  same  situation  holds  for  the  sub¬ 
sonic,  one  atmosphere  exposures,  shown  in  Tables  15  and  16  denoted  by 
Reference  (18)  except  that  30-60  second  exposures  were  performed  for 
the  latter  cases.  Finally,  the  iridium  coated  graphite  and  JTA  tests  in 
Table  14  designated  by  Reference  (17)  were  exposed  for  400-1200  seconds. 
Temperatures  were  measured  with  a  two  color  pyrometer. 

C.  Correlation  of  Results 


Tables  1-16  provide  a  valuable  empirical  guide  to  estimating 
surface  temperatures  from  specific  stream  conditions .  However,  it  is 
desirable  to  provide  a  means  for  comparing  results  contained  in  Tables 
1-6  which  are  related  (i.e. ,  cases  where  a  given  material  is  exposed  to 
two  slightly  different  stream  conditions).  In  other  terms,  it  would  be 
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useful  to  interpolate  and  extrapolate  these  in  order  to  compare  results 
obtained  in  different  facilities  and  to  predict  the  anticipated  tempera¬ 
ture  for  any  given  exposure. 

The  simplest  means  of  performing  such  a  correlation  is  to 
tlciui  wc  the  conditions  for  radiation  equilibrium  at  the  stagnation  point 
of  the  model  on  the  assumption  that  the  energy  lost  by  radiation  is  equal 
to  the  heat  transferred  to  the  model.  On  this  basis, 

CT*T4  =  he<ie'iw[T’  PJ  >  BTU/ft2sec  (1) 

where  <7  =  0.47  x  10'12  BTU/ft2sec  °R^,  t  is  the  total  hemispherical 
emittance,  he  is  the  stagnation  point  heat  transfer  coefficient,  ie  [T,Pe] 
is  the  stagnation  enthalpy,  and  is  the  enthalpy  of  air  at  the  wall  (sur¬ 
face)  of  the  model  and  T°R  is  the  wall  temperature.  Eq.  1  usually  ignores 
reactions  at  the  model  surface  which  give  rise  to  significant  evolution 
or  absorption  of  energy  and  is  as  a  result  but  a  crude  approximation  to 
estimating  the  surface  temperature  which  depends  on  a  knowledge  of  h8. 

If  the  cold  wall  heat  flux,  qcw>  is  defined  as 

'lew  =  Ve  BTU/ft2sec  (2) 

then  Eq.  1  becomes 

T*  =  W1  -  V  [  T*  Pel  /lJ  <3) 


Eq.  3  describes  the  wall  temperature,  T°R  as  an  explicit  function  of  the 
stream  parameters,  ie,  qcw,  and  the  stagnation  pressure  Pe,  as  well 
as  the  enthalpy  of  air  at  the  wall  ^  [  T,  Pe]  and  the  total  hemispherical 
emittance  of  the  surface  c  . 

Numerical  values  of  the  enthalpy  of  air  at  the  wall  are  given 
in  Table  17.  Specification  of  c  ,  qcw,  ie  and  Pe  fixes  the  surface  tem¬ 
perature.  Eq.  3  can  be  solved  numerically  by  using  the  values  of 
U  T,  Pel  shown  in  Table  17.  However,  it  is  convenient  to  represent 
iwl  T,  Pe]  by  analytic  functions  in  order  to  obtain  algebraic  solutions. 
The  following  equations  have  been  employed  for  this  purpose: 


i  [T,  P  ]  *  (T/lOOO)  (100+45  log  P  +(T/l000)  (46-14  log  P  ))  BTU/lb  (4) 

W  C  6  C 


(5) 


i  [  T,  Pe]  =  33.9  (T/1000)2  (2.0  -  0.13  logP)  BTU/lb 

The  numerical  values  of  T,  Pe]  described  by  Eq.  4  and 
Eq.  5  are  compared  with  the  established  values  of  the  enthalpy  of  air  (21), 

(10)  in  Table  17.  In  the  temperature  and  pressure  range  of  interest  (i.  e. . 
2700°  <  T  <  7200°R  and  +?.  ^  leg  P  >  -2),  E4.  4  represents  a  reason¬ 
ably  good  representation  of  i^[  T,  P^|  7  Thus  at  pressures  which  are 
equal  to  or  greater  than  0. 1  atm  the  difference  between  T-  Pe1  and 
Eq.  4  are  less  than  200  BTU/lb.  At  0.01  atm,  larger  differences  are 
noted  at  5400°  and  6300°R.  However,  reference  to  Tables  1-16  shows 
that  under  low  pressure  testing  conditions  ie  is  generally  in  the  5000- 
18,  000  BTU/lb  range.  Under  these  circumstances,  an  error  of  400 
BTU/lb  ini^T,  Pe]  is  not  serious. 

Eq,  5  is  a  poorer  representation  of  i^jT,  Fe]  .  However, 
its  simple  quadratic  form  permits  direct  solution  of  Eq.  3  as  follows: 

0.47«  (T/1000)4  =  qcw  (1-33.9  (T/1000)2  (2.0-0. 13  iogPe)/ie) 

(6) 

hence 

(T/l000)4+33.9q  (2.0-0.13  log  PJ  (0.47«i  )‘ 1  (T/l00O)Z-q  /0.47«  =0  (7) 

cw  ©  ©  cvr 


or 

(T/1000)2  a  -0.5b  +  0.5  (b2-4c) 1/2  (8) 

where 

b  a  33.9qcw  (2.0-0. 13  log  Pfl)  (0.47«  g"1  (9) 

and 

c  3  -*cw/o-47«  (10) 

Eqa  6-10  can  be  employed  to  obtain  a  crude  estimate  of  the 
surface  temperature  at  pressure  equal  to  or  less  than  1  atm.  At  higher 
pressures,  temperatures  computed  to  be  greater  than  54O0°R  will  be 
too  low. 

Eq.  3  has  been  employed  to  compute  for  fixed  values  of 
qCy,  and  T  at  various  stagnation  pressures  and  emlttance  levels.  These 
calculations  were  performed  by  using  the  established  values  of  g[  T,Pe] 
given  by  Reference  (10).  The  results  are  contained  in  Figures  1  to  15 
and  constitute  exact  solutions  to  Eq.  3, 
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Reference  to  Figures  1  to  15  illustrate  the  expected  effect 
of  the  total  omittance  on  surface  temperature.  Thus,  under  stream 
conditions  corresponding  to  a  cnlH  w»Ll  heat  flux  of  2CG  ETE/£t~sec, 
a  stagnation  pressure  of  one  atmosphere  and  a  stagnation  enthalpy  of 
10,  000  BTU/lb;  a  material  having  a  total  emittance  of  0. 2  would  reach 
a  surface  temperature  of  620C°Il  {Figure  6).  under  the  same  condi¬ 
tions,  a  material  having  an  emittance  of  1 . 0  would  reach  a  surface 
temperature  of  only  4400°R  (Figure  10) .  These  curves  also  show  that 
with  increasing  stagnation  enthalpy,  the  surface  temperature  Is  depend¬ 
ent  only  on  heat  flux  and  emittance  as  implied  by  Eq.  3  for  (i^/i aero. 
However,  when  the  heat  transfer  coefficient  is  large,  the  temperature 
depends  upon  the  stagnation  enthalpy.  In  other  terms,  even  under  high 
flux  conditions,  the  surface  temperature  cannot  exceed  the  gas  tempera¬ 
ture.  Note  that  in  all  of  the  proceeding  discussions,  heats  of  chemical 
reaction  at  the  surface  are  ignored.  Naturally,  such  effects  could  re¬ 
sult  in  generation  of  surface  temperatures  which  are  in  excess  of  stream 
temperatures.  The  inflections  present  in  the  le[  T]  curves  at  log  Pe  = 
-2.0  (Figures  11-15)  result  from  the  temperature  dependence  of  the 
enthalpy  as  indicated  in  Table  17. 


Eq.  3  has  been  employed  to  compute  the  surface  tempera¬ 
ture  for  each  of  tho  exposures  shown  in  Tables  1-16  by  employing  the 
representation  of  T,  F]  given  by  Eq.  4.  These  calculations  are 
performed  on  a  computer  using  a  Newton-Raphson  technique  to  obtain 
solutions.  The  resultant  ratio  of  computed  temperature  to  observed 
temperature  is  contained  in  Tables  1-16.  The  ratio,  T(CALC)/T(OBS),  , 
is  plotted  as  a  function  of  cold  wall  heat  flux,  q^,  stagnation  enthalpy, 
i0,  stagnation  pressure,  Pe,  Mach  Number  and  observed  temperature  in 
Figures  16-20.  These  calculations  have  been  performed  using  a  total 
hemispherical  emittance  of  0. 6  for  all  exposures  in  Tables  1  -16  except 
those  in  which  emittance  values  were  measured.  In  the  latter  cases, 
the  reported  values  were  employed. 


An  alternate  calculation  of  the  "Radiation  Equilibrium"  sur¬ 
face  temperature  can  be  performed  which  employs  the  Fay-Riddell 
relation  (22)  to  compute  the  heat  transfer  coefficient,  This  proce¬ 
dure  is  performed  by  setting 


K  a  0.94  ( 


2rB 


du 

« 


1/2 

)  ( 


THT 

6 IV* 


1/2 

> 


(PH) 


1/10 


lbs /it2 


sec 


Letting 


,2RB  ■*“. 

S' 


Q 


for  convenience  and  noting  that 


y-i 


(13) 


Q2  =  * 


gy -  1)  M2  +  2)  [1+  (jy-  1)  M2 4  2)) -I 

(y+  i)  m2  2  -  2  -  J 


(2yM  -  (y  ~  1)) 


for  Mach  Number «  areator  _ i 


Q 


2 


=  9 


(14) 


for  Mach  Number  =  0,  where  y  a  C  /C 
is  the  Mach  Number.  Approximating  v 


1.4  for  an  ideal  gas,  and  M 


((PH)W  /(pp)e)1/10 


(15) 


yields 


he  =  0.94Q1/Z  ((pn)  u„  /2RB)1/Z 


(16) 


24)  the  following: 


Eq.  16  can  be  evaluated  by  estimating  the  viscosity  by  (23, 


H  =  2. 17  x  10"8  T  *^2  lbs  (force)  sec 

iSF 

Setting  P  a  pRT  and  R  =  1724  fi2/seC2  °R  and  noting  that 

<VT«>  *  ^O^ty-IJM2 

withM2  a  uJ/yRT^  yields 


(17) 


(18) 


he  =  0.940 


l/2r 


[1.09x10  8  S6C)  (l/R)l/2M(l-fO.  5(y-  DM2)'1/2]  */2(p  /r  )l/2 

It  xv  c 

(19) 

Also  for  high  Mach  Numbers  (M  >>  1),  Q  a  i.  n  jmdl 

he=  0.99  (6.90  x  10  °)  !/2  (p^/R^)1^2  lb*  (force)  sec/ft3  (20) 
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(21) 


2 

where  Pg  is  in  lbs  (force) /ft  and  Rg  is  in  feet.  Thus 

h  =  2.60  x  10 "5  (,PjRri)1,Z  lbs  (force)  sec/ft3 

6  C  Jj 

(lbs  (force)  =  32,2  lbs  ft/sec  ).  Hence 

he=  8.37  x  10'4  (Pe/RB)1/2  lbs/ft2sec  (22) 

When  the  stagnation  pressure  is  in  atmospheres.  (1  atm  =  2117  lbs  (force)/ft  ) 
he  =  0.0386  (Pg/Rg) 1/2  lbs /ft2 sec  (23) 


Eq.  23  is  a  reduced  Fay-Riddell  relation  applicable  at  high  Mach  Numbers. 

The  general  relation  given  above  in  Eq.  19  is 

he  =  0.0245Q1/Z  [MZ(1  +  0.5  (y  -  1)  M2)"1]  1/4  (Pe/Rg)1^2 

lbs /ft2  sec  (24) 

where  P  is  in  atmosphere  and  Rg  in  feet.  If  the  quantity  0.0245Q  [  1  + 

0.5  (y  -  1)  M^)-!]  1/*  is  approximated  by  0.  0386/(l  +  0. 17M-1),  which  is 
accurate  to  within  5%  for  M  >  0.1,  then  Eq.  24  becomes 

he=  0.0386  (1  +  0.17M'1)'1  (Pe/RB)1/2  lbs/ft2sec  (25) 


Substitution  into  Eq.  1  yields 

<TcT4  =  0.0386  (1  +  0.  17M'1)"1  (24  Pe/D)1/2  (ie-iw[T,Pe])  (26) 


where  D  is  the  diameter  of  the  hemispherical  cap  in  inches. 

Eq.  26  has  also  been  applied  to  compute  the  radiation  equili¬ 
brium  temperature  for  all  of  the  exposures  in  Tables  1  - 1 6  by  employing 
the  description  of  i^  T,  Pe]  afforded  by  Eq.  4.  These  calculations 
were  performed  for  flat  faced  cylinders  by  setting 


D  a  2.  5  Diameter  of  Cylinder 


(27) 
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in  Eq.  26  in  order  to  account  for  the  difference  in  heat  transfer  between 
hemispherical  caps  and  flat  faced  cylinders.  The  stagnation  point  heat 
flux  to  a  flat  faced  cylinder  having  a  diameter  Dc  is  equivalent  to  the 
heat  transfer  to  a  hemisphere  having  an  effective  diameter,  Deff,  equal 
to  fDc.  The  values  of  f  noted  in  the  literature  are  2. 1  (25),  2.5  (12) 
and  2.9  (12).  Tn  the  following  numerical  values  of  i  have  oeen 

reported:  2.0  (26),  3.2  (26),  3.3  (27)  and3.08,  3. 34  and  3. 72  at  Mach 
Number  2.0,  3 .IT and  4.  7^7  respectively  (28).  Thus  it  appears  that 
values  of  2.0  <  f  <  3.72  have  been  employed  for  relating  flat  faced 
cylinders  to  hemispheres.  As  indicated  above,  a  value  of  2.5  is  cur¬ 
rently  being  employed  in  the  present  calculations  (Eq.  27).  This  varia¬ 
tion  may  result  from  individual  facility  characteristics  and  measurement 
techniques.  For  example,  the  value  f  =  2.0  (26)  has  a  heat  transfer 
basis  while  f  a  3.2  (26)  has  a  pressure  gradient  basis.  All  of  the  remain¬ 
ing  values  except  2.T~(25)  were  obtained  experimentally. 

The  results  of  applying  Eq.  26  contained  in  Tables  1-16  under 
the  heading  "Fay-Riddell"  are  displayed  in  Figure  21  which  shows  the 
ratio  T(CALC)/T(OBS)  vs.  the  observed  cold  v/all  heat  flux. 

D.  Discussion  of  Results 

Examination  of  Tables  1-16  and  Figures  16-21  shows  the 
wide  divergence  between  observed  and  calculated  temperatures.  Al¬ 
though  it  is  presently  impossible  to  define  the  causes  of  these  discrepan¬ 
cies,  some  of  the  possibilities  are  worthwhile  noting.  To  begin  with, 
the  exposure  times  for  the  tests  under  consideration  are  variable.  Thus 
the  Avco-MeB2  and  HfB2  +  SiC;  the  Plasmadyne-Zr02  and  JTA;  the 
General  Electric  ATJ,  RVA  and  JTA;  and  the  Cincinnati  Testing  Labora¬ 
tory  JTA  and  Ir  coated  graphite  exposures  were  of  100-1800  seconds 
duration.  The  remaining  tests  were  for  30-60  seconds  with  the  excep¬ 
tion  of  the  CAL  exposures  which  were  limited  to  15  seconds.  Under  these 
conditions,  it  might  be  expected  that  radiation  equilibrium  is  more 
readily  attained  in  the  long  time  exposures  (i.  e. ,  times  greater  than 
100  seconds).  Thus,  the  observed  temperature  in  a  short  time  exposure 
would  be  lower  than  the  calculated  radiation  equilibrium  temperature. 

A  second  major  source  of  error  is  the  measurement  of 
surface  temperature.  This  depends  upon  the  particular  value  of 
spectral  emissivity  employed  in  correcting  the  observed  brightness 
(€  a  1)  to  the  true  temperature.  For  the  G.  E.  and  CTL  results,  ob¬ 
tained  with  a  two  color  pyrometer,  the  error  in  assigning  a  true 
surface  temperature  depends  upon  how  well  grey  body  conditions  are 
approximated.  Although  these  temperature  measurement  errors  can 
be  significant,  it  is  not  likely  that  they  are  the  prime  source  of  the 
present  discrepancies  in  the  T(CALC/OBS)  ratios,  A  related  source 
of  error  is  the  value  of  total  emittance  used  in  computing  the  tempera¬ 
ture  in  Eq.  3,  Aa  indicated  earlier,  measured  values  have  been  employed 
where  available.  Where  no  values  are  available,  a  mean  value  of  0.  6  has 
been  employed.  If  the  true  value  of  total  emittance  is  1.0,  then  the  cal¬ 
culated  temperature  will  be  too  high  by  the  one  fourth  root  of  (1/0.  6)  or 
fourteen  percent.  Such  corrections  could  improve  the  present  results. 
However,  such  changes  would  not  eliminate  the  current  level  of  disagree¬ 
ment  between  observed  and  calculated  temperature. 
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The  relative  sizes  of  the  model  and  the  arc  are  an  additional 
variable  which  has  not  been  considered  in  the  analysis.  Table  13  con¬ 
tains  the  results  obtained  during  the  exposure  of  silicide  coated  re- 

alloy  foil.  In  thccc  the  uiou^oUx  wad 

about  one  inch  and  the  model  was  a  two  inch  square  foil.  The  foil 
radiated  from  both  the  front  and  back  face.  As  a  result,  the  ratio  of 
T ( CALC)  / T ( OB S)  is  much  larger  than  unity.  These  ratios  are  not 
included  in  Figures  16-21. 


Comparison  of  the  T(CALC)/T(OBS)  values  shown  in  Tables 
1-16  for  the  "Cold  Wall"  and  "Fay-Ridden"  Heat  Transfer  Coefficient 
computations  indicate  that  in  some  cases  substantial  differences  exist 
between  these  heat  transfer  coefficients.  Thus,  the  "Cold  WaU"  and 
"Fay-Ridden"  heat  transfer  coefficients  differ  markedly  in  the  data 
generated  at  Avco  and  CTL.  There  is  some  evidence  that  this  may  be 
due  to  turbulent  test  streams  (29).  On  the  other  hand,  the  CAL, 

General  Electric,  Aerospace  Corp. ,  Flasmadyne,  North  American  and 
Grumman  data  show  good  correlation  between  the  "Cold  Wall"  and  "Fay- 
RiddeU"  values.  Reference  to  Figure  21  showB  that  the  ratio  of  T(CALC)/ 
T(OBS)  is  poorer  for  the  "Fay-Riddell"  calculation  than  for  the  "Cold  Wall" 
calculation. 


As  previously  indicated,  the  present  calculations  ignore  heat 
liberated  at  the  surface  due  to  oxide  formation  or  heat  absorbed  by  the 
surface  in  order  to  vaporize  or  melt  the  model  material.  Inclusion  of 
the  former  effects  could  raise  the  present  values  of  T(CALC)/T(OBS) 
while  indue  ion  of  the  latter  would  lower  the  current  T(CALC)/T(OBS) 
ratios.  Although  these  effects  can  be  substantial,  they  do  not  appear  to 
be  the  source  of  the  current  differences.  Inclusion  of  the  heats  of  com¬ 
bustion  would  raise  T(CALC)  and  further  aggravate  the  T(CALC)/T(OBS) 
ratios.  Moreover,  reference  to  Figure  20  shows  that  the  largest  devia¬ 
tions  from  T(CALC)/T(QBS)  =  1,  occur  below  5000°R  where  copious 
vaporization  of  the  materials  under  consideration  is  unlikely. 

The  most  important  (and  unfortunately  the  least  tractable) 
sources  of  error  are  the  reported  values  of  the  stagnation  enthalpy  and 
cold  wall  heat  flux.  Nonuniform  flux  and  enthalpy  conditions  in  addition 
to  gas  radiation  losses  and  variation  in  the  techniques  employed  for  the 
measurement  of  these  quantities  is  probably  the  most  important  single 
source  of  the  deviation  of  T(CALC)/T(OBS)  from  unity  in  Figures  16-21. 
Thus,  comparison  of  theAvco  (half-filled  squares)  and  General  Electric 
(rimmed  squares)  results  for  ATJ  graphite  shows  an  extremely  large 
discrepancy.  TheAvco  data  yield  values  of  T(CALC)/T(OBS)  which  fall 
between  1.10  and  1. 18  based  on  observed  "Cold  Wall"  heat  transfer 
coefficients  while  the  General  Electric  results  for  the  same  material 
range  from  1.  13  to  2.05.  It  should  be  pointed  out  that  the  latter  results 
were  obtained  for  long  time  exposures. 

Apart  from  the  General  Electric  data,  values  of  T(CALC)  / 
(OBS)  which  differ  from  1.0  +  0.2  are  frequent  in  the  Aerospace  and  CAL 
exposures.  As  indicated  earlier,  these  exposures  were  for  short  times 
and  the  high  values  may  be  due  to  the  fact  that  radiation  equilibrium  was 
not  established. 
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E.  Results  Obtained  during  the  Present  Testing  Program 

In  the  course  of  the  present  program,  Loan  /uu  arc 

plaaind.  (.eats  were  performed  (6).  Measurements  of  stagnation  enthalpy, 
stagnation  pressure,  cold  wall'&eat  flux  and  brightness  temperature 
were  performed  in  each  case.  The  latter  were  converted  to  true  tem¬ 
perature  by  using  the  measured  values  of  the  emitta~ce  for  the  candidate 
materials  (5)  under  oxidizing  conditions.  In  addition,  radiated  heat  flux 
was  also  measured  in  these  tests  so  that  normal  emittance,  *  n,  could  be 
deduced  (6).  Table  18  summarizes  the  average  values  of  «  and  the 
T(CALC)/T(OBS)  ratios  for  each  of  the  materials  tested  (6).  The  latter 
ratio  is  computed  on  the  basis  of  Eqs.  3  and  4  corresponding  to  radiation 
equilibria. 


Ideally,  if  radiation  equilibria  were  the  dominant  factors  and 
all  measurements  were  accurate,  these  ratios  should  be  unity.  Although 
there  are  departures,  it  is  satisfying  to  note  that  the  differences  are 
small  compared  to  those  obtained  by  considering  the  results  of  other 
studies  (i.  e..  Figures  16-21  and  Tables  1-16)*  Reference  to  Table  18 
shows  that  ratios  of  T(CALC)/T(OBS)  are  lower  for  cases  where  melting 
is  observed  than  for  cases  where  a  solid  oxide  (or  coating)  is  present. 
Moreover,  Table  18  shows  that  large  values  of  T(CALC)/T(OBS)  are 
characteristic  for  some  of  the  materials.  The  occurrence  of  ratios  which 
are  larger  than  unity  implies  resistance  to  energy  absorption  by  the 
material.  Thus,  exposure  of  H£B£  +  SiC(A-4)  and  HfC  +  C(C-ll)  to 
identical  stream  conditions  (i.  e. ,  stagnation  pressure,  enthalpy  and  cold 
wall  heat  flux)  would  result  in  an  1 1%  lower  surface  temperature  than  that 
reached  by  HfC  +  C(C-ll).  This  conclusion  would  apply  if  stream  condi¬ 
tions  were  not  sufficient  to  produce  melting  of  HfBg  +  SiC(A-4).  At  lower 
levels,  KT-SiC(E-14),  WSi2/W(G-18)  and  Sn-AZ /Ta-10W(G-19).  which 
exhibit  T(CALC) /T(OBS)  ratios  of  1.43,  1.56  and  1.41,  respectively 
demonstrate  similar  resistance  to  energy  transfer.  Although  the  origin 
of  this  resistance  is  not  clear  at  present,  it  iB  probably  due  to  blocking 
effects  caused  by  evolution  of  gaseous  oxides.  These  observations  sug¬ 
gest  a  method  of  ranking  the  behavior  of  the  refractory  materials  which 
differs  from  the  customary  recession  vs.  temperature  curves.  In 
Section  IV,  an  alternative  method  of  presentation  which  compares  reces¬ 
sion  rate  as  a  function  of  heat  flux  and  enthalpy  for  the  candidate  materials 
is  considered.  This  method  does  not  require  a  knowledge  of  the  spectral 
or  the  normal  emittance  and  integrates  the  blocking  effects  characteristic 
of  each  material. 
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HI.  SIMPLIFIED  MODEL  FOR  CALCULATING  THE  OXIDATION 

bENavior  oP  Graphite  under  Nigh  velocity  aIR  STOW 

CONDITIONS' - 

A.  Introduction,  and  Summary 

In  the  course  of  the  present  study,  the  oxidation  of  graphites 
in  air  has  been  investigated  experimentally  over  a  range  of  conditions  (4- 
f>)  between  2500°?-  and  6500°R,  at  velocities  between  i  ft/ sec  and  Mach  3.  2. 
Additional  studies  (12),  (13)  have  been  performed.  The  succeeding  volume 
of  this  series  (30)  presents  a  complete  discussion  of  the  surface  reaction 
problem  encountered  in  the  oxidation  of  graphite.  This  discussion  con¬ 
siders  the  coupling  of  mass  transport  through  the  boundary  layer  with 
reactions  at  the  surface  in  detail  (30).  By  contrast,  the  discussion  pre¬ 
sented  below  employs  simplified  models  (7,  31  -33)  which  provide  an 
explicit  means  for  computing  the  rate  of  grapEIte  recession  as  a  function 
of  density,  surface  temperature,  gas  velocity,  stagnation  pressure  and 
sample  radius.  The  results  are  compared  with  observations  covering  a 
range  of^ensity  between  80  and  115  lbs/ft^,  temperatures  between  2500° 
and  6500°R,  velocities  between  1  ft/sec  and  Mach  8.0,  stagnation  pres¬ 
sures  between  0.007  and  1. 0  atm  and  nose  radii  between  0.005  and  0.  07  ft. 

B.  Derivation  of  the  Simplified  Model  for  Graphite  Oxidation 

Comparison  of  the  recession  rates  observed  for  graphite  in 
the  present  study  at  air  flow  rates  between  1  ft/sec  and  250  ft/sec  (Figure 
37  and  pp  23-29  of  reference  5)  demonstrate  that  most  (if  not  all)  of  the 
graphite  oxidation  data  previously  determined  (34,  35)  is  controlled  by  a 
supply  limit.  Thus,  for  example,  Eqs.  28  (GuIEransen  et  al.  (34)) 
corresponds  to  a  recession  rate  of 

rti  a  1.86  x  10  F0  e  gm/cm  sec  ( 

°2 


where  PA  is  in  torr  and  T  is  in  *K,  and 
u2 

*  •>  aa  „  m-^0'  32  -3600/RT  , 

S  =  2.44  x  10  e  mils/  min 

2 

for  graphite  with  a  density  p  =  1.80  gms/cm3  or  112  lbs/ft3.  For  air,  Po2 
=  160  ,  thus  Eq.  28  yields  a  rate  of  0.  0007mils/secat  3000°R.  This 
compares  with  values  near  0. 1  mils/ sec  at  a  flow  rate  of  1  ft/ sec  and  1 
mil/ sec  at  250  ft/sec  observed  in  the  current  study  (5,  Figure  37).  The 
oxidation  rates  reported  by  Okada  and  Ikeqawa  (35)  between  1800°R  and  5400°R. 
under  slow  flow  conditions  at  a  pressure  of  0.21  atm  02  yield  similar  com¬ 
parisons.  Values  corresponding  to  0.02,  0.07  and  0. 10  mils/sec  are 
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reported  (35)  at  temperatures  of  2000°R,  3000°R  and  4000°R,  respectively. 
As  indicate?  (5),  transport  of  sufficient  quantities  of  oxygen  to  the  reacting 
surface  and  alcnowledge  of  the  oxygen  partial  pressure  at  the  reacting  sur- 

face  arps  rfifliiirftH  in  nrHpr  fn  malrA  ™Cd"i,ngful,  ThcCIC- 

tical  studies  along  these  lines  are  presented  in  reference  (30)  in  order  to 
calculate  concentration  gradients  present  in  the  gaseous  boundary  layer 
adjacent  to  the  reaction,  surface.  However,  since  the  details  of  this  treat¬ 
ment  are  quite  complex  it  is  instructive  to  consider  the  approximate  treat¬ 
ment  described  earlier  (Eq.  15  p.  242,  Reference  (7)). 

o  The  correct  picture  for  the  oxidation  of  graphite  above  800°C 

(1472°F)  is  one  of  continually  increasing  rate  with  both  temperature  and 
oxygen  pressure  according  to  the  Arrhenius  relationship: 


ih  =  kPz  e-E/RT 


From  800°  to  2065°C  (1472°-3750°F)  z  =0.32-0.38  and  E  a  3600-4200 
cal/mol  based  on  results  of  Gulbransen  (34)  and  this  study  on  Speer  710 
and  RVA  graphites  (5).  Although  the  activation  energies  observed  in  the 
present  study  (5)  are  comparable  to  those  reported  by  Gulbransen  et  al. , 
the  rates  are  much  higher  in  the  present  investigation.  Thus,  for  P02  = 
150  torr,  at  T  =  1700°F,  Eq.  28  yields  a  rate  of  about  0.04  mils/minute. 
The  present  study  indicates  60  mils/minute  at  250  fps.  Eq.  28 
converts  to 


(29) 


m 


=  3. 16  x  10'5  P°’  32  e"3600/RT  lbs/ft2 


sec 


(30) 


where  is  in  atmospheres  and  T  is  in  °K,  and 
u2 

S  =  3.39  *10-3p"-32„-3600/RT  mU./ 


sec 


(31) 


for  graphite  with  a  density  p=1.80  gms/cm3  or  112  lbs/ft^.  These  reac¬ 
tion  kinetic  equations  are  altered  (7)  to  reflect  diffusion  control  (31)  yielding 


m 


=  k(CQ  M/Mq  )ZPZ  (l-m/mD)ze'E/RT  lbs/ft2 
2,  e  2 


sec 


(32) 


where  Pe  is  tl  e  stagnation  pressure  and  C02  e  i®  the  maos  fraction  of 
oxygen  at  the  edge  of  the  boundary  layer.  ' 

If  one  considers  the  reaction 


C  +  l/202  +  2N2-*  CO  +  2N2 


(33) 


as  dominant,  then  M.  =  28.5,  Mo?  =  32  and  C q2  a  0.21.  In  addition  (7, 
32)  6 
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xxx ^  --  u.uuo  iba/tt  sec  (34) 

Thus,  Eq#  32  describes  the  oxidation  of  graphite  in  terms 
of  an  Arrhenius  term  (ke“E/RT)  and  a  pressure  correction  term  (1-ih/ihj*)) 
which  is  based  on  a  diffusion  limit  mQ  given  by  Eq.  34.  In  other  terms, 
the  (l-ih/mjj)  coefficient  can  be  considered  as  a  pressure  correction  which 
relates  the  oxygen  concentration  at  the  edge  of  the  boundary  layer  to  that  at 
the  graphite  surface.  Divergent  values  have  been  reported  for  the  Arrhenius 
constants  k  and  E(5).  Part  of  the  difficulty  undoubtedly  arises  from  the  fact 
that  CO?  is  the  dominant  product  gas  of  graphite  oxidation  at  low  temperature 
while  CO  dominates  at  high  temperature  (36).  In  contrast  to  the  values  k  = 

3.  16  x  10*5  and  E  =  3600  implied  by  Eq.  5 (T,  Scala  and  Gilbert  (32)  have  pro¬ 
posed  k  =  6. 729  x  108  and  E  =  44, 000  for  "fast  kinetics"  and  k  ="4. 473  x  104 
and  E  -  42, 300  for  "slow  kinetics".  The  current  results  indicate  that  lower 
values  of  the  activation  energy  may  be  more  appropriate  (5). 

As  indicated  below  the  results  of  the  present  Btudy  (5,  30)  have 
been  examined  in  order  to  obtain  the  most  appropriate  valueB  of  k  andTT,  To 
begin  with,  Eq.  34  was  re-examined  in  order  to  allow  for  air  flow  rate  ef¬ 
fects  at  s iib sonic  velocities  since  the  diffusion  limited  rate  for  carbon  re¬ 
moval,  mg3,  was  originally  approximated  by  (7) 

*KD  *  A*Pe/RB>‘/2cK  (35> 

where  Gr  is  the  mass  fraction  of  carbon  in  the  oxidation  reaction.  Since 
the  average  molecular  weight  in  Eq.  33  is  25.8  (with  Mo?  =  32)  then  = 

1/7.  The  coefficient  "A"  was  estimated  by  analogy  with  the  Fay- Riddell 
relation  (22,  27) 

q  =  0.042  (P^Rg) 1/2  (i^)  (36) 

where  the  enthalpy  difference  (ie-iw)  is  the  analog  •'i  Cv  in  the  mass  loss 
relation.  Thus,  with  A  =  0.042,  m^  is  defined  by  I*),  34.  However,  at 
low  velocities  a  Mach  Number  correction  (of  Eq.  25)  is  required  and  the 
result  yields  (30) 

mD  =  0.0072  (1  +  0.17M"1)"1 

with 

*  .  0.74P‘,3.-10'730'RT 

u2 

for  the  reaction  rate,  then 


(pe/RB)1/Z  lbfl/ft2sec  (37) 


lbs/ft2sec  (38) 
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ifa  =(0. 170)l/3  [1-mtRg)1^2  /0.0072  {1  +  0. i7M-1)-1  F^1/Z]  1/3 

-10  /*»•*•  .  z 

t,w-c.  74  c  - *  lbs /ft  sec  (39) 

where  R-,  is  the  body  radius  in  feet.  Since 
£> 

S  (mils/sec)  =  m  (Ibs/ft^aec)  12000/p  (lbs /ft3)  (40) 


Hence 


S  =  4920p_1pi/3  exp(-9720/T) 
6 


,l/2«-l/2  1/3 


0.0U6p§R"£P^1,i  J 

[  i - — rrrr- 1 

(1  +  0.17M  1 


mils/  sec  (41) 

for  the  carbon  recession  rate  in  air,  where  p  is  the  density  of  graphite  in 
lbs/ft3  and  T  is  in  °R.  The  activation  energy  and  pre-exponential  factor 
in  Eq.  38  were  evaluated  by  employing  the  form  of  Eq.  41  and  the  data  in 
Table  19  in  order  to  obtain  the  beat  fit. 

The  experimental  results  shown  in  Table  19  include  results 
for  ATJ,  RVA(B-5),  PT0178(B-9)  and  Poco  Graphite  (B-10).  The  calcula¬ 
tions  refer  to  fiat  faced  cylinders  (Rg  =  2.  5  Rc)  exposed  during  the  present 
study  by  R.  A.  Tanzilli  (13),  and  by  Metzger  et  al.  (12).  Exposures  per¬ 
formed  in  the  current  atu3y  are  designated  by  run  numbers  (61.  Tanzilli' s 
exposures  of  PT0178(B-9)  axe  designated  by  B9-Gel-18  in  Table  19.  The 
exposures  of  ATJ  graphite  by  Metzger  et  al.  are  denoted  by  MED.  In 
addition  to  these  arc  plasma  exposures,  calculations  were  performed  for 
sexveral  high  velocity  CG/HW  tests  conducted  on  RVA(B-5)  (5).  Exposure 
B5-L2,  B5-L3  and  B5-.L4  were  taken  to  represent  the  results  contained  in 
Figure  45  of  (5)  at  200  ft/ sec.  Examination  of  the  samples  after  exposure 
indicates  a  nose  radius  equal  to  one  sixteenth  of  an  inch  as  shown  in  Figure 
44  of  (5),  hence  Rjj  =  62  mils  or  0. 0052  ft.  The  data  and  calculations 
cover  temperatures  between  2300°  and  659Q°R.  Thus,  it  is  possible  that 
some  of  the  high  temperature  points  (T  >  6000°R)  reflect  recession  via 
diffusion  limited  vaporization  where  the  observed  values  would  be  ex¬ 
pected  to  exceed  Eq.  41  since  vaporization  becomes  a  significant  factor 
in  the  recession  rate  above  6000“R  at  one  atmosphere.  The  present  data 
cover  the  pressure  range  between  0.007  and  1  atmosphere  and  represent 
RVA(B-5),  PT0179(B-9)  and  Poco(B-lO)  graphite.  Moreover,  the  results 
include  data  generated  in  the  HG/CW  Avco  Model  500  and  ROVERS  facilities 
(3)  and  the  General  Electric  Tandem  Gerdien  Arc  as  well  as  the  Lockheed 
M/S  Co.,  CG/HW  facility  (2). 
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Ia  addition,  it  should  be  pointed  out  that  an  independent 
experimental  arc  plasma  study  by  Sallis  et  al.  (3?)  between  2500°  and 
5000°P.  Crt  stagflation  hohveen  0. 1  and  9  atmosDheres  yields 

a  pre -exponential  of  3.  0  and  an  activation  energy  of  15,  500  cal/ mole 
as  compared  with  the  current  results  of  0.74  and  10,730  cal/mole  shown 
in  Eq.  28.  In  addition,  Sallis  et  al.  suggest  that  the  exponent  of  the 
oxygen  pressure  is  0. 4  rather  than  the  value  of  0.333  indicated  by  Eq.  38. 


Eq.  41  has  been  employed  to  generate  Figures  22-26, 
which  compare  with  the  observed  and  computed  results  of  graphite 
oxidation  at  low  velocities.  Figure  22  shows  excellent  agreement  for 
RVA(B-5)  at  3310°R  (Figure  37  of  reference  5),  while  Figure  23 
shows  qualitative  agreement  (see  Figure  45  61  reference  5).  While  it 
is  apparent  that  Eq.  41  provides  an  excellent  descriptionfor  the  oxida¬ 
tion  kinetics  of  graphite  over  a  very  wide  range  of  conditions  shown  in 
Table  19,  it  is  of  interest  to  consider  the  results  of  a  different  kind  of 
experiment  relative  to  the  predictions  of  Eq.  41. 


Blyholder  and  Eyring  (38)  have  measured  the  oxidation  of 
graphite  in  flowing  oxygen  at  800°K.  "The  oxygen  pressure  employed  was 
26  microns  of  mercury  while  the  flow  rate  was  1000  cm/  sec.  The  sam¬ 
ples  employed  in  this  experiment  were  hollow  cylinders  of  carbon  having 
a  density  of  1. 3  gms/emy  (81  lbs/ ft^)  which  were  1/4  inch  in  diameter 
with  a  one  millimeter  wall  thickness.  The  samples  were  cut  in  half 
parallel  to  their  longitudinal  axes  prior  to  exposure.  Under  these  condi¬ 
tions,  the  oxidation  rate  of  carbon  (as  measured  by  formation  of  CO) 
corresponded  to  16  x  10^  molecules  of  CO/cm^sec.  This  corresponds 
to  3. 2  x  10““  gms  of  carbon/ cm2 sec.  In  order  to  apply  Eq.  41  to  compute 
the  rate  of  carbon  oxidation  for  comparison  with  this  result  is  is  necessary 
to  estimate  the  value  of  the  term  within  the  braces  in  Eq.  41  firBt.  Since 
the  flow  rate  was  1000  cm/sec  or  about  33  ft/ Bee,  M  is  approximately  0.03. 
The  corresponding  value  of  Fe  is  equal  to  3.4  x  10"5  atm.  (26p  of  Hg)  mul¬ 
tiplied  by  (1.00/0.21)  or  160  x  10““  atm.  Estimating  Rj$  ^  2.5  Re  loads 
to  Rq  a  0. 026  ft.  Thus,  Eq.  41  becomes 


S  a  3.85  x  10'3  (1-80&)1/3  mils/sec  ( 

or  §  =  3.4  x  10"3  mils/sec  a  8.  6  x  10"^  cm/ sec  a  11.2  x  10"^  gms/cm^sec. 
This  value  is  in  good  agreement  with  the  observed  value  of  3,2  x  10““  gm/ 
cm2sec  in  view  of  the  estimates  required  for  Rq  and  M.  Moreover,  this 
experiment  waa  performed  far  from  the  range  of  conditions  employed  to 
fix  the  pre -exponential  and  activation  energy  in  Eq.  38. 

C.  Comparison  of  the  Scala-Gilbert  and  John-Schick  Models 
for  Graphite  Ablation 

John  and  Schick  (33)  have  developed  a  theory  for  describing 
the  diffusion  controlled  ablation'ol  graphite  which  defines  the  linear  reces¬ 
sion  rate,  §(mils/sec),  of  graphite  as  follows: 


IV, 


INFLUENCE  OF  FLUX -ENTHALPY  AND  ALTITUDE- VELOCITY 

VARIABLES  ON  THE  TOJtJE^lOTTOF'REFRA'CrOKY  MATERIALS 


Figures  1-8  of  reference(6)  describe  the  30  minute  recession  or 
oxidation  depths  observed  for  the  candidate  matmriaim  s«  HG/CW  arc 
plasma  tests  and  in  CG/HW  furnace  tests  as  a  function  of  surface  tempera¬ 
ture.  In  the  former  case,  the  surface  temperature  is  a  tesult  of  the 
interaction  of  the  material  with  the  stream.  Thus,  while  comparison 
of  a  given  material  in  the  CG/HW  and  HG/CW  cases  based  upon  results 
at  a  given  temperature  is  quite  legitimate,  evaluation  of  various  materials 
in  HG/CW  tests  on  solely  a  temperature  basis  is  not  complete.  As  indi¬ 
cated  earlier  (6)  factors  such  as  emittance,  oxidation  products  and  surface 
characteristics  can  lead  to  situations  where  identical  stream  conditions 
produce  a  variety  of  surface  temperatures  on  different  materials  oven 
after  long  exposure  times,  hi  order  to  consider  an  alternative  method  for 
comparison  of  the  oxidation  characteristics  and  to  relate  the  HG/CW  testa 
to  flight  parameters,  an  additional  description  can  be  employed. 

Figure  27  shows  the  stagnation  pressure  as  a  function  of  altitude 
and  velocity.  In  addition,  stagnation  heat  flux  to  a  one  inch  sphere  described 
as  a  function  of  altitude  and  velocity  Is  displayed.  In  essence,  these  curves 
provide  a  means  for  relating  the  HG/CW  arc  plasma  tests  to  flight  trajectory 
conditions.  The  relations  are  presented  for  the  case  of  a  one  inch  radius 
sphere.  Since  the  heat  flux  is  proportional  to  R"1  /2,  the  heat  flux  to  a  4 
inch  sphere  would  be  one  half  of  the  values  shown  on  the  ordinate  of  Figure 
27,  while  the  heat  flux  to  a  1/4  inch  sphere  would  bo  twice  thoao  shown  in 
Figure  27.  The  curves  relating  stagnation  pressure,  PQ,  altitude,  A,  and 
velocity,  V,  in  Figure  27  are  the  results  of  complex  equations  (39-41). 
However,  between  50  Kilo  feet  and  2S0  Kilo  feet,  these  curves  can  He 
repr e so ntod  simply  by  Hfcq.  !>0  as: 

P0  «  V2  (1  +  (A/216)4)  10 "A/54  atm  (50) 


where  P.  is  the  stagnation  pressure  behind  a  normal  shock  in  atmospheres, 

V  is  the  velocity  In  kilo  feet/sc  e  and  A  is  the  altitude  in  kilo  foot.  The  stag¬ 
nation  enthalpy,  ie,  is  approximately 

ie  «  20V2  BTU/lb  (51) 


Eqs.  50  and  51  when  coupled  with  the  Fay-RiddeU  relation  for  a  hemisphere 

q  »  0.042  i^Pg/Rjj)1^2  BTU/ftZeec  (52) 

where  Rg  is  the  body  radius  iu  feet  yields 

q  a  2.9V3  (1  +  (A/216) 4)1/2Rb"1/210"A/108  (53) 

whoro  Rg  is  the  body  mdius  in  inches. 
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Thus,  Figure  27,  or  the  approximations  afforded  by  Eqs.  50-53 
permit  direct  conversion  from  altitude-velocity  apace  to  flux-enthalpy 
space  for  a  hemisphere.  As  indicated  above,  the  behavior  of  th«  re¬ 
fractory  materials  of  interest  uuJer  KC/C;w  conditions  simulating  high 
velocity  atmosphere  flight  is  presented  in  Figures  1-8  of  reference  6. 

The  latter  show  material  recession  as  a  function  of  stagnation  pressure 
and  surface  temperature.  In  this  case,  the  surface  temperature  results 
from  a  combination  of  q,  ie  and  material  reactions  with  the  stream.  The 
relations  between  q,  ie,  Pe  and  surface  temperature  (for  radiation  equili¬ 
brium)  is  illustrated  in  Figures  1  to  15.  Calculations  of  this  type  have 
been  peri  nned  for  eachHG/CW  exposure  (6).  Although  reasonable  agree¬ 
ment  of  computed  T  [  i@,q,  P  J  with  observe?  surface  temperature  has  been 
encountered,  some  systematic  differences  have  been  noted  as  shown  in 
Table  18.  Specifically,  SiC  and  SiC  bearing  composites  reach  lower 
surface  temperatures  for  given  ie,  q  and  Pe  conditions  than  do  the  other 
materials  considered.  Similar  behavior  is  noted  for  WSig/W  and  SnAl/Ta- 
10W  under  conditions  where  the  coatings  do  not  fail. 

An  alternative  method  for  describing  the  performance  of  these  re¬ 
fractory  materials  is  illustrated  in  Figures  28-36.  These  figures  show 
heat  flux  as  the  ordinate  and  stagnation  enthalpy  as  the  abscissa.  In 
addition,  velocity  as  related  to  stagnation  enthalpy  by  Eiq.  51  is  shown 
as  the  abscissa.  Moreover,  the  relationship  for  a  one  inch  radius  sphere 
at  150  kilo  feet  ie  the  ordinate.  Thus,  for  a  velocity  of  16  kft/sec  at  an 
altitude  of  150  kft,  the  stagnation  enthalpy  would  be  approximately  5100 
BTU/lb  and  the  stagnation  heat  flux  to  a  one  inch  radius  sphere  would  be 
about  600  BTU/ft2sec.  If  the  body  radius  were  4",  the  heat  flux  would  be 
300  BTU/ftzsec  (located  by  dropping  down  to  300  on  the  inner  ordinate 
acale).  Conversely,  if  the  body  radius  were  1/4",  the  heat  flux  would  be 
1200  BTU/ftrsec  (located  by  moving  up  to  1200  on  the  inner  ordinate  scale). 
Thus,  the  double  set  of  ordinate  and  abscissa  scales  permit  direct  trans¬ 
lation  of  velocity  and  body  radius  to  flux  altitudes,  Eq.  53  has  been  em¬ 
ployed  to  construct  the  inserted  curve  shown  on  each  graph  which  shows 
the  ratio  of  qf  Altitude]  /  q[  15C  kft]  as  a  function  of  altitude.  This  ratio  is 
0.27,  0.48,  1.00,  2.46  and  6.85  at  250,  200,  150,  100  and  50  kft, 
respectively.  Thus,  at  a  velocity  of  16  kft/sec  and  250  kft  altitude,  the 
heat  flux  to  a  one  inch  sphere  would  be  600  x  0. 27  or  162  BTU/ft2sec. 

Under  these  conditions  a  4"  radius  would  experience  a  heat  flux  of  81 
BTU/ft2sec  while  a  1/4"  radius  sphere  would  be  exposed  to  a  heat  flux 
of  324  BTtl/ft2sec.  Similarly  at  an  altitude  of  50  kft,  the  heat  flux  to  a 
one  inch  radius  sphere  would  be  2.46  x  600  =  1476  BTU/ft^eec  and  the 
flux  to  a  4"  and  1/4"  sphere  would  be  738  and  2952  BTU/f^sec,  respectively. 

Thus,  Figures  28-36  show  heat  flux  and  enthalpy  for  any  velocity, 
altitude  and  body  radius.  Figure  28  shows  the  recession  rates  observed 
for  hafnium  diboride  at  stagnation  pressures  of  one  atmosphere  (circles) 
and  0.01 -0.1  atmospheres  (squares).  Recession  rates  which  are  less  than 
0. 1  mils/sec,  between  0. 1  and  1  mil/ bbc  and  more  than  one  mil/  sec  are 
indicated  by  open,  half- filled  and  filled  poiuts,  respectively.  Recession 
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rates  are  plotted  at  flux  and  enthalpy  co-  ordinates  for  exposures  given 
in  reference  6.  At  present,  sufficient  data  are  not  available  to  con¬ 
struct  boundaries  representing  constant  recession  levels  over  the  entire 
flux- enthalpy  space.  To  bridge  the  gap,  1.0  mil/sec  boundaries  are 
approximated  at  each  pressure  by  temperature  levels  obtained  from 
Figures  1-15.  Figure  36  summarizes  all  of  the  results  for  a  1  mil/  sec 
boundary  at  one  atmosphere. 

hypsrbclic  curves  £cir  stU  of  m o.tori^ls  exespt  ^rMphitcs 
and  tungsten,  define  the  flux-enthalpy  (or  velocity- altitude-body  radius) 
conditions  where  the  recession  rate  exceeds  1  mil/ sec.  Flux-enthalpy 
conditions  below  and  to  the  left  of  these  boundaries  result  in  recession 
rates  which  are  less  than  1  mil/  sec.  Flux-enthalpy  conditions  above  and 
to  the  right  of  these  boundaries  yield  recession  rates  in  excess  of  1  mil/ 
sec.  In  the  case  of  graphites  and  tungsten,  the  linear  boundaries  (based 
on  diffusion  limits)  are  shown  for  0.  5,  1.0,  2.0  and  4.0  mils/ sec.  Flux- 
enthalpy  conditions  lying  below  and  to  the  right  of  these  boundaries  will 
result  in  lower  recession  rates,  while  those  lying  above  and  to  the  left 
will  result  in  higher  rates,.  The  computed  rates  for  graphite  and  tungsten 
are  calculated  on  the  basis  of  Reference  31  and  Eq.  45. 

Although  Figure  36  is  based  ona  limited  number  of  tests,  it  provides  a 
clear  indication  of  the  superiority  of  SiC  and  SiC  composites.  Figure  1 
of  reference  6  suggests  that  ZrB2(A-3)  exhibits  recession  rates  below 
1  mil /  sec  at  surface  temperatures  up  to  5000°F.  Figure  5  of  Reference 
6  shows  that  KT-SiC(E-14)  exhibits  rates  below  1  mil/ sec  or  less  below 
?000°F.  However,  Table  18  as  well  as  Figure  36  indicate  that  the  flux- 
enthalpy  conditions  which  produce  5000°F  surface  temperature  (and  a  1 
mil/ sec  recession  rate)  for  ZrB2(A-3)  will  yield  a  surface  temperature  of 
only  4000°F  (and  a  comparable  recession  rate)  for  KT-SiC(E-14).  Of 
course,  relative  mechanical  properties,  thermal  shock  resistance,  density 
and  other  factors  may  impose  additional  criteria  for  comparison.  Never¬ 
theless,  Figure  36  presents  a  direct  ranking  of  oxidation  behavior  as  a 
function  of  flight  conditions  for  extended  periods  of  time  at  the  stagnation 
point.  The  position  of  these  curves  could  vary  with  stagnation  pressure. 
There  are  several  other  inversions  of  rank  relative  to  the  CG/HW  and  HG/ 
GW  tests  using  temperature  as  a  base.  Figures  42  and  45  of  reference  6 
show  that  HfB2  +  SlG(A-4)  and  (A-?)  exhibit  1  mil/ sec  recession  rates  aF 
lower  temperatures  than  HfC  +  C(C-ll).  Nevertheless,  for  given  stream 
conditions,  the  latter  reaches  surface  temperatures  which  are  20%  higher 
than  the  former  (see  Table  18)  and  ranks  lower  in  Figure  36.  Similarly, 
Si02  +  W(H-22)  and  Si/RVC(B-8)  which  degrade  rapidly  at  4000°F  and 
3100QF  respectively  in  GG/HW  furnace  tests  rank  high  on  the  basis  of 
Figure  36.  The  relative  behavior  of  Si/RVG(B-8)  and  Ir/C(I-24)  is  also 
illustrated  graphically  in  Figure  36,  Reference  to  Figures  5  and  8  of 
reference  6  shows  that  Si/RVC(B-8)  provides  protection  for  graphic  at 
surface  temperatures  up  to  37O0°F  in  the  HG/CW  arc  plasma  tests  while 
Ir/C(I-24)  is  protective  up  to  surface  temperatures  of  4200°F.  However, 
since  the  latter  las  an  emittance  of  0. 3  and  T(CALC)/T(OBS)  =1.16 
(Table  18)  while  the  former  has  an  emittance  of  0.70  and  T(CALG)/T(OBS)  = 
1,43,  Si/RVA(B*8)  exhibits  a  much  greater  resistance  to  heat  flux  and 
enthalpy  than  does  Ir/C(I-24).  Addition  of  Hf02  to  the  Ir/C(I-24)  coating 
system  improves  the  performance  of  this  material  by  increasing  its 
emittance  ((>). 
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V 


.  UTILIZATION  OF  THE  FLUX-ENTHALPY  VS.  ALTITUDE- 
VELOClTY 

SPECIE  IEEHPRXJ  ECTC5IE5 - 


In  order  to  i1l>i«?trato  the  wum  by  which  Figure  36  can  be  employed 
to  predict  the  behavior  of  a  candidate  refractory  material  for  a  specific 
trajectory  it  is  worthwhile  to  consider  the  Air  Force  Flight  Dynamics  Laboratory 
FD1-7MC  lifting  reentry  vehicle's  maximum  cross  range  characteristics. 

This  vehicle  is  designed  for  a  Lift/Drag  ratio  between  2.5  and  3.0. 

Table  20  provides  the  altitude  and  velocity  as  a  function  of  time. 

These  data  were  employed  to  calculate  stagnation  enthalpy,  pressure  and 
heat  flux  based  on  Eqs,  50-53  for  a  3"  nose  radius.  The  results  are  shown 
in  Table  20  and  in  Figures  37-40.  These  figures  indicate  that  HfB2+SiC 
(A-4)  and  (A-7),  HfB?(A-2),  ZrB2+SiC(A-8)  and  HfC+C(C-ll)  could  survive 
the  entire  trajectory  in  this  configuration.  In  addition,  ZrB2(A»3)  and 
ZrB2+SiC+C(A- 10)  might  also  survive.  Comparison  of  the  flux-enthalpy 
values  for  this  trajectory  with  the  results  of  arc  plasma  tests  shows  that 
HfB2(A-2),  ZrB2(A-3),  HfBz+SiC(A-4)  and  (A-7),  ZrB2+SiC(A-8),  and 
ZrBj+SiC+CJA-lO)  survived  tests  which  are  equivalent  tp  the  FDL-7MC 
with  very  small  recessions^  Due  to  the  low  temperature  oxidation  of 
HfC+C(C-l  1),  this  material  might  be  limited  for  reuse  (6).  The  borides 
and  boride  composites  would  not  suffer  from  this  limitation.  A  number  of 
long-time  cyclic  exposures  of  diboride  composites  have  been  performed  (6) 
in  the  Model  500  and  ROVERS  facilities  to  evaluate  reuse  capabilities  for  “* 
trajectories  typified  by  FDL-7MC.  The  results  provide  a  striking  illus¬ 
tration  of  the  reuse  capability  of  these  materials  for  lifting  reentry  applica¬ 
tions  . 


Sample  i+20%SiC(A-7)-28R  was  exposed  for  thirteen  cycles  at 

0.07  atm  stagnation  pressure,  a  stagnation  enthalpy  of  10,200  BTU/lb  and 
a  cold  wall  heat  flux  of  495  BTU/ft2aec.  Each  cycle  was  about  1800  seconds 
long  with  a  total  exposure  time  of  22,500  seconds.  The  surface  temperature 
increased  from  one  cycle  to  the  next  starting  at  3500°R  and  holding  near 
5300°R  for  cycles  5  through  13.  Total  material  recession  was  15  mils  after 
this  extremely  long  exposure.  Sample  ZrB2,  l+20%SiC(A-8)-15M  was  exposed 
for  four  cycles  at  1.0  atm  stagnation  pressure,  a  stagnation  enthalpy  of  5000 
BTU/lb  and  a  cold  wall  heat  flux  of  380  BTU/ft^sec,  Each  cycle  wag  1800 
seconds  long,  total  exposure  time  was  7200  seconds.  The  surface  tempera¬ 
tures  were  near  5000°R.  Total  material  recession  was  26  mils.  Finally, 
sample  ZrBZ+SiC+C(A-10)-26R  was  exposed  at  0.  236  atmospheres  stagna¬ 
tion.  pressure,  a  stagnation  enthalpy  of  7700  BTU/lb  and  a  cold  wall  heat  flux 
of  455  BTU/ft^sec.  This  test  covered  eleven  cycles  of  approximately  1800 
seconds  duration  for  a  total  exposure  time  of  18,900  seconds.  Surface 
temperature  held  near  5100°R  after  the  first  cycle.  Total  material  reces¬ 
sion  was  83  mils. 


These  results  illustrate  the  reuse  capability  of  boride  composites  for 
lifting  reentry  application,  since  they  exceed  the  range  of  conditions  and 
FDL-7MC.  This  capability  is  unrivaled  by  any  other  materials  system. 
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VI.  CALCULATION  OF  THE  FLUX-ENTHALPY  BOUNDARIES  FOR 


Figures  28-40  display  the  recession  rates  observed  in  arc  plasma 
tests  z.  function  of  Leal  uu&  and  stagnation  entnaipy.  This  representa¬ 
tion  indicates  the  locus  of  heat  flux,  q,  and  stagnation  enthalpies,  i.  e., 
which  define  the  region  where  recession  rates  exceed  1  mil/  sec  for 
refractory  materials.  This  value  was  arbitrarily  chosen  in  order  to 
illustrate  a  means  by  which  the  flux-enthalpy  representation  can  be  em¬ 
ployed.  The  1  mil/ sec  boundaries  can  be  identified  by  collecting  suf¬ 
ficient  data  to  cover  the  q-ie  space  completely.  Since  sufficient  data 
points  are  not  available  to  do  this  experimentally,  the  procedure  em¬ 
ployed  in  Figures  28-35  is  to  associate  the  1  mil/  sec  boundary  for 
condensed  oxide  forming  refractory  materials  with  a  specific  temperature. 

This  procedure  is  quite  arbitrary. 

An  alternative  method  is  to  consider  the  1  mil/  sec  rate  as  being 
characteristic  of  melting.  Under  these  conditions,  the  convective  heat 
flux  can  be  considered  as  a  source  of  the  radiative  losses  and  the  heat  of 
melting.  This  heat  balance  is  defined  by  Eq.  54  as  follows: 

he(ie-iw[T,P])  =  <T«  (T/lOOO^-MlZ.OOOrW^  BTU/ftZsec  (54) 


where  the  heat  transfer  coefficient,  h  ,  is 

© 

“e-  lew'*,  <55> 

The’  enthalpy  of  air  at  the  wall,  [T,  P]  ,  1b  given  by  the  following  expres¬ 
sion  from  Eq.  (4), 


ijT,P]  =  (T/1000)  (100  +  451ogP  +  (T/1000)  (46-141ogP))  BTU/lb  (56) 

and  q,~»  and  i.  are  the  cold  wall  heat  flux  and  stagnation  enthalpy,  respec¬ 
tively! 

In  Eq.  54,  8*0.47  BTU/ft2sec°R4  is  Boltzmann's  constant,  T°R 
is  the  melting  point  of  the  refractory  material,  p  is  the  density  in  lbs/ fir, 
c  is  the  total  normal  omittance,  AHj  is  the  latent  heat  of  melting  in  BTU/lb 
and  &  is  the  recession  rate  in  mils/  sec.  Since  there  is  not  available 
experimental  data  for  AHf,  estimates  have  been  made  as  indicated  in 
Table  21.  Thus,  Eq.  54  Is  similar  to  the  radiation  equilibrium  surface 
temperature  calculation  except  that  an  additional  heat  loss  (12,000"*pAH£S) 
has  been  included.  Rearrangement  yields: 
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mils/ sec 


g  =  (q/V  ^e'VtT.P])  -  0,  (T/IOOO)4 
(12,  000) -1  pAHf 

where  fte’eo"  ££? hSS  “■***<« 

“3*??.  to  melt  the  material  ~a',  r^eTte^.Tien  r"e“'‘ 

The  latent  -Tj*  ‘ AHi  aa^  p.  for  4116  candidate  materials, 

since  virtually  no*  mea^urementsYre ^aiTablef  FigL^l^fah'ow^18 

EfyJdE  SUSS  20-35*  “  agreement  £«h  tie  r^«^” 


VH,  CALCULATION  OF  TEMPERATURE  DISTRIBUTION  THROUGH  A 

OTMF05ITE  CYLINDER  UNDER  STEADY  STATE  CONDITIONS - 

gCEOWTOCTm ;  5IDET05SE5~VIA  'R7tDIAT  ION - 

The  calculation  of  surface  temperature  from  stream  conditions  based 
on  radiation  equilibrium  described  by  Equations  (1-4)  in  Section  VII  has  been 
performed  for  all  teste.  The  results  wLR.Ii  die  given  in  Table  io  show 
that  reasonable  agreement  can  be  obtained  in  most  cases  to  within  10  or  15%. 
In  most  cases,  the  calculated  temperature  is  too  high.  At  present,  the  re¬ 
sults  indicate  that  the  over  estimated  values  are  due  to  a  reduction  in  heat 
transfer  coefficient  due  to  vaporizing  oxidation  products.  This  conclusion 
is  based  on  the  fact  that  materials  containing  SiC,  SiO?  and  other  vapor 
products  (i.  e.  Sn-Af/Ta-10W(G-19)  yield  higher  values  of  T(CALC)/T{OBS) 
as  seen  in  Table  18.  The  following  calculation  deals  with  the  problem  of 
side  losses  via  radiation  and  axial  temperature  gradients  in  the  arc  plasma 
teat  cylinders. 

Representation  of  the  steady  state  temperature  distribution  through  a 
right  circular  cylinder  of  radius  R,  length  L,  and  a  coating  which  has  a  thick¬ 
ness  I  as  shown  below  can  be  defined  on  the  basis  of  a  convective  heat  input 
and  radiation  losses  from  the  front  and  sides  according  to  the  model  shown 
above.  In  this  description,  the  thermal  conductivity  of  the  coating  and  sub¬ 
strate  are  k.  and  k  ,  respectively. 

1  s 


No  Heot  Loss 
Via  Conduction 


The  total  heat  balance  requires  that 

(q/ie)  (ie-iw(Tf,  PJ)  =  ff«  T*  <1  +  (2L/R)  (Fp  +  Fg))  BTU/ft2sec  (58) 

where  q  and  i  are  the  cold  wall  heat  flux  and  stagnation  enthalpy  respectively, 
and  the  enthalpy  at  the  wall,  i  ,  at.  a  temperature  T,  and  a  stagnation  pressure 
P  is  approximately  (from  Eq.  (4)), 

i  [Tf,F  ]  =  (Tf/1000)  (100  +  45  log  Pfi  +  (Tf/l000)(46-14  log  Pfl))  BTU/lb 

(59) 
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Evaluation  of  the  temperature  distributions  along  the  length  of  the  cylinder 
are  performed  on  the  basis  of  the  following  assumptions, 

(a)  the  gradient  through  the  coating  is  linear  for  0  <  x  <  I 

(b)  the  temperature  distribution  through  the  substrate  is  quadratic 
fnr  T  <  y  <■  T_ 

(c)  heat  losses  at  x  =  L  are  negligible 

(d)  there  are  no  radial  temperature  gradients. 


With  these  assumptions,  the  quantities  Fp  and  Fg  contained  on  the  right 
side  of  Eq.(58)  are  defined  as  follows 

Si 

FF  =  J  (T/Tf)4  dx  feO) 

o 


where  x  =  x/L,  xj  s  I/L  and  T  is  the  temperature  at  any  value  of  x.  Similarly, 


(«s/«)  (T/Tf)4  dx 


(61) 


where  <§  is  the.  emittance  of  the  substrate  and  <  is  the  emittance  of  the  coating. 
If  the  ratio  T/T|  is  defined  as  T,then 

T  =  1  -  (x/x^  for  0  <  x  <  Xj  (62) 

where  is  equal  to  T/T£  at  x  »  I,  and 

T  =  (l-xj)"2  [  Ti  +  (Tt  -  (2  -  xj)]  for  x,<  x<  1  (63) 


where  Tt  =  Tt/Tr  is  the  ratio  of  the  back  face  temperature  Tt  to  the  front  face 
temperature  Tj.  Matching  of  the  conductive  fluxes  at  the  interface  requires  that 

Si 


Tj/Tf  =  T £  .  Ux^Rd-ify]  ■1[2^11?t+k(l-3ei)Csx1(l-Il)  JVdx]  (*1) 


o 
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where  k  =  kp  /kg  is  the  ratio  of  thermal  conductivities  of  the  coating  and  sub¬ 
strate  and  Cg  is  defined  by  Eq.i(65)  as  follows 


Cq  =  2«<TT?L2/Rk 


(65) 


{62 )3,nd  (^4/pcriiiii  cvaludUuu  u£  jc  since 


f  4dx  =  ^  (1  +  x(T .  -  l)xi“1)4  dx 


Ff  =  S.  [  1  +  2<T.  -  1)  +  2<T.-1)2  +  (T.-l)3  +  (1/5)  (T.-l)4], 
Similarly,  defining 

a3  =  tfi  -  ftV2-Si>l  (i"5!)”2  anda4  =  (T^T.Jaj1  (l-x^’2 

yields 

(«/*s)Fg  =  J  a4  [  1  +  a4x(2-x)]  4  dx 


1.0 

5 

h 


4 

=  times 


a4  *  S' +  xi  I  ’  T  *?  +  *?  "  T1  *i  +  'T  *i  “  TT  ] > 


4x4[^^.Rf+^Si.2]) 


+  a2(i^  +  2 xJ[-4-  x2  +  3^-4]) 


f  a4  <  *T  +  4*i  [  T  *i  "  XJ ) 


+  1  -  x, 


(66) 


(67) 


(68) 


(69) 
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Eq.  (43)can  be  employed  to  define  the  temperature  gradient  in  the  cylinder 
at  the  front  face 


I  dT  \  _  /,  jr  »—  -1 

!dr!*  =  ° - 


(70) 


Since  all  of  the  energy  entering  the  cylinder  at  its  front  face  must  be 
radiated  away,  then 


/  dT  \ 

( 1=~  ^  * 
dx  x  =  0 


1-T 


-  Cs*'1  S  I4  ®  =  csr*  (FF+FS> 


(71) 


or 


(Ff  +  Fs)  »  (l-T.Jk/xjCg  (72) 

'When  =  I  =  0,  i.  e, ,  no  coating  is  present,  Eq,  (62)  yields 

-(“?->  =  2(1-TJ  b  CgFg  (73) 

dx  x  a  0  1  a  “ 


or 


Fs  =  W-TJ/Cg 


(74) 


The  limiting  case  where  I  »  L  «  0,  indicating  no  radiation  losses  from  the  sides, 
reduces  Eq,(58)  to 


(q/ie)  (ie  -  iw  [Tf,  Pe])  *  CT«T£4  (75) 

for  radiation  equilibrium  from  the  front  face.  This  result  is  Identical  to 
Equations  1  -4  of  Section  II. 
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Solution  of  Eq.  (58)  for  the  general  case  where  I  and  L<  are  not  equal 
to  zero  requires  location  by  iteration  of  the  proper  value  of  T,  which  satis¬ 
fies  Eqs.  (58),  (64),  (73),  (68)  and  (72).  This  procedure  is  started  by 
employing  the  solution  of  Eq.  (58)  (based  on  1  =  L  =  0)  for  T,  as  a  first 
guess.  Once  this  first  guess  is  available,  C  is  defined  by  Eq.  (65).  These 
initial  values  of  T..  and  C  are  emoloved  in  conjunction  with  Eqs .  (64).  (66) 
and  (68)  to  search^for  a  vllue  of  7,  satisfying  (72).  The  iteration  procedure 
is  begun  by  substituting  the  resultant  value  of  (Fp  +  F„)  into  Eq.  (58)  and  solv¬ 
ing  for  the  resultant  error.  The  procedure  is  rcpcatca  by  perturbing  the  ini¬ 
tial  Tp  repeating  the  solution  and  obtaining  the  resultant  error  for  the  second 
choice.  Examination  of  these  values  of  T^  and  the  resultant  errors  ,  using  a 
Newton- Raphson  method,  permits  the  iteration  process  to  proceed  to  conver¬ 
gence  . 


When  1  =  0  the  procedure  is  the  same  except  that  Eq.  (72)  is  replaced 
by  Eq.  (64).  In  this  case  Fp  is  zero.  A  computer  program  has  been  developed 
to  perform  these  calculations. 

Sample  calculations  are  shown  in  Table  22  for  arc  plasma  tests  ZrB? 
(A-3)-2MC  which  exhibited  a  front  face  temperature  of  4930°R  and  an  internal 
temperature  of  3400  R  at  a  distance  of  100  mils  from  the  front  face  (6).  The 
effects  of  the  internal  pyrometer  hole  was  considered  theoretically  and  found 
to  be  negligible.  Consequently,  the  present  calculations  are  based  on  the  total 
length  of  429  mils.  The  results  illustrate  the  large  temperature  gradients 
through  the  sample  indicated  in  reference  6. 

In  order  to  gain  some  insight  into  the  effects  of  length,  L,,  radius,  R, 
coating  thickness,  1,  and  the  thermophysical  properties,  «  ,  «  g,  k_.  and  In  on 
the  temperature  gradients,  it  is  instructive  to  reconsider  Equatioft  (58).  in 
Equation  58,  the  calculation  of  T.  reduces  to  the  simple  radiation  equilibrium 
case  given  by  Equation  75  and  Equations  1-4  of  Section  11  when  L  is  equal  to 
zero  or  when  R  is  infinitely  large  compared  to  L.  When  this  is  not  the  case, 
the  temperature  gradients  are  controlled  by  the  values  of  Fp  and  Fg  in 
Equation  58.  These  quantities  are  in  turn  dependent  on  the  related  values  of 
T,  and  T  defined  in  Equations  62-64.  The  relation  between  T,  and  T,  (Equa¬ 
tion  64)  specifically  involves  the  thermophysical  properties  through  Cg  given 
by  Equation  65.  The  explicit  correspondence  between  the  radiation  parameters 
F„  and  Fg  and  the  thermophysical  parameters  T. ,  T.  and  Cg  are  given  by 
Equations0  72  or  74,  Figure  45  shows  the  variation  of  TT  ana  T,  specified  by 
Equation  64  for  a  representative  case  where  ^-./kg  »  0.10.  Figure  46  illus¬ 
trates  a  portion  of  the  coating  contribution  (F Jj  to  the  radiation  term  in  Equa¬ 
tion  58  and  its  dependence  on  T. ,  The  total  radiation  parameter  {Fp  +  Fg) 
varies  with  T  as  shown  in  Figure  47,  This  dependence  is  illustrated  for  the 
case  where  i/L  =0,10  and  kp/k_  a  0.  10  on  the  basis  of  Equations  66  and  69. 
Solutions  for  T.  and  (F„,  +  Fg)  are  given  by  the  indicated  intersections  using 
Equation  72.  These  calculations  are  carried  out  by  the  computer  program 
which  then  solves  Equation  58.  In  the  special  case  where  1=0  (no  coating 
present),  Fp  vanishes  and  Fg  as  given  by  Equation  69  is  solely  dependent  on 
T  as  shown  in  the  figure  based  on  Equation  74. 
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The  effect  of  the  thermophysical  properties  and  the  radiation  losses 
on  the  surface  temperature  gradient  within  the  material  is  evident  from  Equa¬ 
tion  71.  For  small  Cq  the  temperature  gradient  is  small,  as  may  arise  for 
low  surface  temperature,  (T,),  physically  thin  cylinders  (L-*  0),  relatively 
thin  cylinders  (L/R  -*  0),  or  Targe  thermal  conductivity  (k  )  aft  of  any  coating. 
For  large  values  of  those  parameters,  C„,  and  thus  the  temperature  gradient, 
tend  to  be  large.  Relatively  small  thermal  conductivity  of  the  coating  (i.e. 

T;  -  k_/kg  a, Via  11)  alnu  Icado  iu  larger  fnermui  gradients,  iaistly,  the  nor¬ 
malized  radiation  parameters,  F_,  and  influence  the  gradient  by  partially 
compensating  for  the  CQ  effect.  Figure  47  shows  that  small  (F„  +  Fn)  is 
associated  with  large  and  vice  versa.  Physically,  this  implies  larger 
gradients  at  the  face  lead  to  lesser  radiation  losses  from  the  aft  parts  of  the 
cylinder  due  to  the  lesser  temperatures  then  present. 

Tables  23-28  summarize  the  results  obtained  by  comparing  the  observed 
internal  temperatures  (6)  with  calculated  values  based  on  Eqs,  58-75  for  ZrB^+ 
SiC(A-8) ,  ZrB2(A-3),  l!iB2+SiC(A-7),  RVA(B-5),  ZrB2+SiC-t-C(A-10),  WSi2/W 
(G-18)  and  Hf-Ta-Mo(I-23).  These  tables  contain  the  measured  front  face 
temperature  Tp,  the  observed  temperature  at  a  distance,  d  mils  from  the  front 
face,  T(d),  and  the  cold  wall  heat  flux,  q,  the  stagnation  enthalpy,  ie,  and  the 
stagnation  pressure,  Pe.  In  addition,  these  tables  show  the  radius,  R,  length, 

L,  and  oxide  coating  thickness,  I,  The  latter  was  equated  to  the  conversion 
depth  for  the  oxide  formers  (6).  For  WSi2/W,  I  was  equated  to  the  WSi2  coat¬ 
ing  thickness  and  1=0  for  RVX(B-5)  graphite  which  ablates  without  coating 
formation,  Values  of  the  cmittancc,  «g,  taken  from  Table  18  as  well  as  suit¬ 
able  values  of  the  thermal  conductivities  of  the  thermal  conductivity  of  the 
coating  kjr  and  the  substrate  ks  are  also  shown  In  Tables  Z3-Z8. 

The  computed  results  are  displayed  In  terms  of  the  ratio  of  calculated 
front  face  temperature  to  observed  front  face  temperature  Tf  (CALC) /Tf  (OBS) 
and  the  ratio  of  computed  In  depth  temperature  T^(CALC)  to  computed  front 
face  temperature  T>  '(CALC).  If  the  agreement  Is  exact  (e.g. ,  Hf-Ta-Mo(I-23)- 
43R  In  Table  28),  the  ratio  of  Tf  (CALC)/Tf  (OBS)  would  be  1.00.  In  the  exam¬ 
ple  Tf  (CALC)  is  4440°R  vs,  4530°R  a  T f  (OBS),  Similarly  the  measured 
temperature  at  120  mils  Is  3560°R  vs.  Td  (CALC) o3380°R.  In  this  case,  the 
observed  gradient  is  960°R  while  the  calculated  gradient  is  1060°R  In  120  mils, 

All  of  the  runs  shown  In  Tables  23-28  were  performed  on  flat  faced 
cylinders  except  those  designated  by  a  suffix  H  (hemisphere)  or  S  (cylindrical 
shroud  with  a  200  mil  wall).  Photographs  of  these  models  havo  been  presented 
(6).  The  shrouds  and  hemispherical  caps  did  not  alter  the  gradients  observed 
for  flat  faced  cylinders.  Thus  all  of  the  calculations  were  based  on  flat  faced 
cylinders  ignoring  the  hemispherical  caps  and  the  shrouds.  Reference  to  Tables 
23-28  indicate  relatively  good  agreement  between  calculation  and  observation, 
in  view  of  the  simple  model  employed  and  the  complexities  of  the  experiments. 

The  largest  deviations  occur  at  low  surface  temperatures  (i.e. ,  Tf  < 
3300°R)  for  the  materials  which  form  SIO2  as  an  oxidation  product.  Thus,  in 
cases  where  samples  of  HfB2+S>iC( A-7),  ZrB2+SiC(A-8),  ZrB2+SiC+C(  A- 10) 
or  WSl2/W(G-18)  were  exposed  with  shrouds  or  as  large  diameter  hemispheres 
Tf  (CALC)  is  considerably  larger  than  Tf  (OBS).  However,  this  difference  is 
smaller  than  obtained  when  Tf  is  computed  on  the  banis  of  rudiation  equilibrium 
(6)  (i.e.,  Eq.  75).  The  causo  of  this  behavior  is  presently  unknown  (6).  Refer¬ 
ence  to  Tables  23-28  shows  that  the  calculated  and  observed  ratios  ofT^/T- 
are  in  general  agreement. 
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D  is  the  hemispherical  cap  diameter  of  the  model.  As  asterisk  indicates  a  flat  face  cylinder  having 
a  geometrical  diameter  of  D  inches. 


SUMMARY  OF  PUBLISHED  FLUX- TEMPERATURE  DATA  FOR  REFRACTORY  MATERIALS 
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D  is  the  hemispherical  cap  diameter  of  the  model.  An  asterisk  indicates  a  flat  face  cylinder  ha' ing 
a  geometrical  diameter  of  D  inches. 
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is  the  hemispherical  cap  diameter  of  the  model.  An  asterisk  indicates  a  flat  face  cylinder  having 
geometrical  diameter  of  D  inches. 


SUMMARY  OF  PUBLISHED  FLUX- TEMPERATURE  DATA  FOR  REFRACTORY  MATERIALS 
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D  is  the  hemispherical  cap  diameter  of  the  model.  An  asterisk  indicates  a  flat  face  cy Linde  r  having 
a  geometrical  diameter  of  D  inches. 
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a  geometrical  diameter  of  D  inches. 
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D  is  the  hemispherical  cap  diameter  of  the  model.  An  asterisk  indicates  a  flat  face  cylinder  having 
a  geometrical  diameter  of  D  inches. 


SUMMARY  OF  PUBLISHED  FLUX- TEMPERATURE  DATA  FOR  REFRACTORY  MATERIAL 
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TABLE  16 


TABLE  IV 


TEMPERATURE  AND  PRESSURE  DEPENDENCE  OF  THE 
ENTHALPY  OF  AIR  (BTU/LB)* 


T 

Log  P 

°n. 

I-?. 

+  1 

0 

-1 

-2 

— 

— 

"  ’  1  "" 

1  1 

— 1  - 

■ 

2700 

704 

704 

704 

704 

704 

707 

707 

709 

709 

709 

(644) 

(625) 

(605) 

(585) 

(566) 

(430) 

(462) 

(494) 

(526) 

(558) 

3600 

968 

968 

968 

971 

980 

980 

980 

981 

983 

991 

(917) 

(937) 

(956) 

(975) 

(995) 

(764) 

(821) 

(879) 

(936) 

(993) 

4500 

1238 

1244 

1263 

1318 

1485 

1274 

1279 

1295 

1346 

1497 

(1220) 

(1300) 

(1382) 

(1462) 

(1543) 

(1194) 

(1284) 

(1373) 

(1462) 

(1551) 

5400 

1530 

1574 

1707 

2055 

2641 

1596 

1634 

1750 

2063 

2630 

1  „ 

(1551) 

(1716) 

(1881) 

(2046) 

(2211) 

(1720) 

(1848) 

(1977) 

(2106) 

(2234) 

6300 

1867 

2053 

2485 

3060 

3293 

1967 

2115 

2489 

3047 

3341 

(1911) 

(2183) 

(2456) 

(2727) 

(3000) 

(2341) 

(2516) 

(269D 

(2866) 

(3041) 

7200 

2321 

2746 

3329 

3596 

3719 

2416 

2759 

3301 

3631 

3789 

(2301) 

(2703) 

(3104) 

(3506) 

(3908) 

(3058) 

(3286) 

(3515) 

(3743) 

(3972) 

8100 

2890 

3495 

3879 

4073 

4523 

2943 

3447 

3877 

4131 

4585 

(2720) 

(3274) 

(3828) 

(4382) 

(4936) 

(3870) 

(4159) 

(4448) 

(4738) 

(5027) 

9000 

3536 

4086 

4382 

4913 

6468 

3515 

4037 

4412 

4964 

6469 

(3168) 

(3897) 

(4626) 

(5354) 

(6084) 

(4778) 

(5135) 

(5492) 

(5847) 

(6206) 

♦  The  sequence  of  valueB  in,  Reference  (21)  Reference  ( 10) 
(Equation  (4)  )  and( Equation  (5)). 
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TABLE  19  (CONT) 


COMPARISON  OF  CALCULATED  AND  OBSERVED  RECESSION 
RATES  FOR  RVA(B-5),  PT01781B-9).  pnrn/n.  1  n\  A  TvT  A  rn  T  /—  *  T^-r-r-r - 

'  '  -  -  *  —  x_4  ivrfciiii  1 

(FLAT  FACE  CYLINDERS) 


Observed 
Sample  Recession 
No-  Rate  _T 

°R 


_ £_  Density 

(atm)  (lbs/ft  j 


B5-23M  0.878  3725  1.000  112 
B5-24M  0.992  4105  1.000  112 
B5-25M  0.522  3420  1.000  112 

B5-26M  0.269  3035  1.000  112 

B5-27M  0.430  2995  1.000  112 

B5-31M  0.538  3285  1.000  112 

B5-32M  0.422  3475  1.000  112 

B5-16R  1.463  5855  0.218  112 

B5-28R  0.018  2165  0.005  112 

B5-29R  0.086  2780  0.008  112 

B5-30R  0.244  3465  0.011  112 


Effective  S(CALC) 
w  Radius  Calculated 

Mach  Diameter  R  Recession 
Initial/Final  B  Rate 


(mils) 

0.  10  502/350 

0.  13  502/330 

0.15  503/400 

0.15  503/380 

0.15  502/340 

0.10  501/350 

0.  10  502/405 

3.20  739/620 

3.20  503/500 

3.20  504/490 

3.20  502/470 


ft.  (mils/sec) 

0.0443  1.275 

0.0433  1.525 

0. 0449  1.420 

0.0460  1.192 

0.0438  1.175 

0.0443  1.174 

0.0472  1.199 

0.0707  1.266 

0.0522  0.074 

0.0517  0.188 

0.0505  0.298 


S(OBS) 

5(CALC1 


0.689 

0.650 

0.368 

0.226 

0.366 

0.458 

0.352 

1.156 

0.243 

0.457 

0.819 
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TABLE  20 


CALCULATION  OF  FLUX-ENTHALPY  CONDITIONS  CORRESPONDING 
TO  THE  TRAJECTORY  FOR  THE  FDL-7MC  HYPERSONIC  RE-ENTRY 
VEHICLE  UNDER  MAXIMUM  CROSS  RANGE  CONDITIONS 

(Rb  =  3") 


Time 

Altitude 

(sec) 

(kft) 

200 

260 

400 

20  0 

600 

J  90 

800 

180 

1000 

180 

1200 

180 

1400 

170 

1600 

160 

1800 

150 

2000 

140 

2200 

130 

2400 

120 

2600 

110 

2800 

100 

3000 

80 

3200 

60 

Velocity 
kft/  e  ec 

26 

25 

23 

22 

20 

18 

16 

14 

13 

12 

10 

8 

7 

5 

3 

2 


BTU/lb 

BTU/ft^sec 

(atm) 

13500 

210 

0.03 

12500 

480 

0.21 

10600 

450 

0.25 

9700 

480 

0.34 

8000 

360 

0  28 

6500 

260 

0.23 

5100 

210 

0.24 

39  00 

170 

0.28 

3400 

170 

0.36 

2900 

160 

0.43 

2000 

110 

0.45 

1300 

70 

0.42 

1000 

60 

0.48 

500 

30 

0.38 

200 

10 

0.30 

100 

5 

0.31 

W/S  *  53  lbs/ft2 

W/CLS  a  312  lbs/ft^ 
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TABLE  21 


SUMMARY  OF  DATA  USED  IN  THEORETICAL  CALCULATIONS  OF 
RECESSION  RATES  4S  A  FUNCTION  OF  FLUX-ENTHALPY  CONDITIONS 


Material 

TV 

£ 

AHf 

Density 

(  8ms  ) 

'  g.  at. 

ASf 

Code 

(BTlJ/lb) 

(cal/ g. 

H£B2>1(A-2) 

6570 

0,40 

490 

625 

66.7 

5.0 

ZrB2(A-3) 

HfB2+20v/oSiC. 

6335 

0.37 

840 

350 

37.6 

5.0 

(A-4) 

5700 

0.48 

495 

585 

57.5 

5.0 

Boride  Z(A-5) 

5300 

0.55 

705 

355 

37.5 

5.0 

H£B2  j(A-6) 
HfB2+20v/oSiC 

6570 

0.50 

490 

665 

66.7 

5.0 

(A-7) 

5700 

0.48 

480 

565 

59.0 

5.0 

ZrB~+20v/ oSiC 

(A-8) 

5300 

0.55 

740 

340 

35.8 

5.0 

H£B,+35v/oSiC 

(A-9) 

5500 

0.59 

560 

485 

49.3 

5.0 

ZrB2+SiC+C 

(A-10) 

5500 

0.62 

1150 

280 

28.7 

6.0 

HfC+C(C-l  1) 

6210 

0.55 

730 

565 

59.5 

7.0 

ZrC+C(C-12) 

5725 

0.45 

1050 

340 

38.0 

7.0 

JTA(D-13) 

5000 

0.57 

2000 

190 

17.5 

7.0 

JT0992(F-15) 

5800 

0.60 

1400 

290 

28.8 

7.0 

JT0981{F-16) 

5170 

0.51 

1880 

195 

19.2 

7.0 

WSi2/W(G-18) 

6570 

0.28 

70 

1200 

184.0 

2.0 

Sn-Ai/ T  a-W  (G-l  9) 

5890 

0.46 

65 

1055 

181.3 

2.0 

H£-Ta-Mo(I-23) 

4320 

0.54 

50 

840 

179.0 

2.0 

Ir/C(I-24) 

4950 

0.30 

50 

1400 

192.0 

2.0 

*  Estimated  values, 

**  Calculated  value,  AH^  =TAS^. 
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TABLE  22 


CALCULATION  OF  TEMPERATURE  GRADIENTS  THROUGH 
ARC  PLASMA  TEST  SAMPLE  ZrB2(A-3) -2MC 


Input  Values 

q  =  365  BTU/ft2sec 

ie  =  3230  BTU/Ib 

Pe  =  1.060  atm 

« S  =  «  =0.4? 

R  a  0.020  ft.  (246  nails) 

L  a  0.035  ft.  (429  mils) 

I  =  0. 0011  ft.  (14  mils) 

kF  =  0.0001  BTU/ft  sec°R  for  oxide 

kg  a  0. 0120  BTU/ft  sec°R  for  boride 


Observed  Values 

Tf  a  4930°R 
T  a  3400°R  at  x  =  101  mils 


X 

Computed  Temperature 

X 

Computed  Temperature 

mils 

°R 

mils 

°R 

0.0 

4633 

231 

3325 

21.0 

3424 

273 

3314 

42.0 

3411 

315 

3306 

63.0 

3398 

357 

3300 

105.0 

3376 

399 

3297 

126.0 

3366 

420 

3297 

168.0 

3347 

COMPARISON  OF  OBSERVED  AND  CALCULATED 
INTERNAL  TEMPERATURE  RESULTS 
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COMPARISON  OF  OBSERVED  AND  CALCULATED 
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Figure  1.  Surface  Temperature-Heat  Flux  Relations  As  A  Function  ef  Stagnation  Pressure,  Enthalpy 
And  Surface  Emittance. 
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Figure  8.  Surface  Temperature -Heat  Flux  Relations  As  A  Function,  of  Stagnation  Pressure,  Enthalpy 
And  Surface  Emittance. 
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Figure  14.  Surface  Temperature-Heat  Flux  Relations  As  A  Function  of  Stagnation  P: 
And  Surface  Emittance. 
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Figure  16.  Ratio  of  Calculated  to  Observed  Surface  Temperatures  Based  on 

the  Cold  Wall  Heat  Transfer  Coefficient  and  Radiation  Equilibrium 
Vs.  Cold  Wall  Heat  Flux. 
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Figure  17.  Ratio  of  Calculated  to  Observed  Surface  Temperatures  Based  on 

the  Cold  Wall  Heat  Transfer  Coefficient  and  Radiation  Equilibrium 
Vs.  Stagnation  Enthalpy. 
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Figure  18.  Ratio  of  Calculated  to  Observed  Surface  Temperatures  Based  on. 

the  Cold  Wall  Heat  Transfer  Ct efficient  and  Radiation  Equilibrium 
Vs.  Stagnation  Pressure, 
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Figure  19.  Ratio  o£  Calculated  to  Observed  Surface  Temperatures  Baaed  on 

the  Cold  Wall  Heat  Transfer  Coefficient  and  Radiation  Equilibrium 
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Figure  20,  Ratio  of  Calculated  to  Observed  Surface  Temperatures  Based  on 

the  Cold  Wall  Heat  Transfer  Coefficient  and  Radiation  Equilibrium 
Vs.  Observed  Temperature. 
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Figure  21.  Ratio  of  Calculated  to  Observed  Surface  Temperatures  Based  on 

the  Fay-Riddell  Heat  Transfer  Coefficient  and  Radiation  Equilibrium 
Vs.  Cold  Wall  Heat  Flux. 
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Velocity ,  lift /see  At  150  kft  For  A  One  Inch  Radius  Sphere 


5. 


Figure  36.  Ideation  of  1  mil/sec  Recession  Rato  Boundary  us  a  Function  oi 
Heat  Flux  and  Enthalpy  for  Refractory  Material*  at  P  «=  l  utm. 

Various  Recession  Levels  are  Shown  for  RVA  Graphite  and 
Tungoton. 
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V'l  >,iire  37.  Receeiilon  Rate  of  Hafnium  Diboride  as  a  Function  of  Heat  Flux 
and  Stagnation  Enthalpy,  Compared  with  FDI..-7MC-MCR 
Trajectory. 
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Figure  38.  Recession  Rate  of  Zirconium  Diboride  as  a  Function  of  Heat  Flux 
and  Stagnation  Enthalpy.  Compared  with  FDL-7MC-MCR 
Trajectories. 
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Velocity,  kft/sec  At  150  kft  For  A  One  Inch  Radius  Sphere 
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Figure  39.  Recession  Rate  of  HfB^  +  SiC  Composites  as  a  Function  of  Heat 
Flux  and  Stagnation  Enthalpy.  Compared  with  FDL.-7MC-MCR 
Trajectory. 
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Velocity .  kft  /tec  At  150  kft  For  A  One  Inch  Radius  Sphere 
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Figure  41.  Computed  Boundaries  for  Melting  Recession  of  HfB^  ^  i(A-2), 

10b 


Velocity  t  kft/sec  At  150  kft  For  A  One  Inch  Radius  Sphere 
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Figure42.  Computed  Boundaries  for  Melting  Recession  of  ZrB2(A~3). 
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Figure  43i  Computed  Boundaries  for  Meltiug  Recession  of  HfC+C  (C-ll). 
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Velocity,  kft/sec  At  ISOkft  For  A  One  Inch  Radius  Sphere 
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Figure  45.  Variation  of  T,  with  Tt  for  a  Number  of 

According  to  Equation  64  when  I/L  »0  arv 

kp/kg  a  0,  10. 
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Figure  46.  The  Effect  of  Interface  Temperature,  T.,  on  the  Level 
of  the  Radiation  Contribution,  Fp,  According  to 
Equation  66, 
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Figure  47.  The  Dependence  o£  the  Total  Radiation  Parameter 

(SV+Fe)  on  for  various  value*  of  Cg  when  I/L  =  0.  10 
ana  kp7kg  «  0,  10  derived  from  Equation*  66  and  69. 
Intersections  Define  Solution*  on  the  Bails  of  Equation 


Figure  48.  The  Effect  of  Cg  on  Tt  and  the  Corresponding 
Radiation  Parameter  Fg  for  the  case  where 
I  a  0  (l.  e.  =  ?£)  Based  on  Equation  74. 
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